Transport and magnetic properties of GdBaCo 2 05 +:r single crystals: A cobalt oxide 
with square-lattice C0O2 planes over a wide range of electron and hole doping 
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Single crystals of the layered perovskite GdBaCo2 0s+2, (GBCO) have been grown by the floating- 
zone method, and their transport, magnetic, and structural properties have been studied in detail 
over a wide range of oxygen contents, < x < 0.77. The obtained data are used to establish a rich 
phase diagram centered at the "parent" compound GdBaCo2 0s.5 - an insulator with Co ions in the 
3+ state. An attractive feature of GdBaC^Os+z is that it allows a precise and continuous doping of 
C0O2 planes with either electrons or holes, spanning a wide range from the charge-ordered insulator 
at 50% electron doping (x = 0) to the undoped band insulator (x — 0.5), and further towards 
the heavily hole-doped metallic state. This continuous doping is clearly manifested in the behavior 
of thermoelectric power which exhibits a spectacular divergence with approaching x=0.5, where it 
reaches large absolute values (±800 /xV/K) and abruptly changes its sign. At low temperatures, the 
homogeneous distribution of doped carriers in GBCO becomes unstable, as is often the case with 
strongly correlated systems, and both the magnetic and transport properties point to an intriguing 
nanoscopic phase separation into two insulating phases (for electron-doped region) or an insulating 
and a metallic phases (for hole-doped region). We also find that throughout the composition range 
the magnetic behavior in GBCO is governed by a delicate balance between ferromagnetic (FM) 
and antiferromagnetic (AF) interactions, which can be easily affected by temperature, doping, or 
magnetic field, bringing about FM-AF transitions and a giant magnetoresistance (MR) phenomenon. 
What distinguishes GBCO from the colossal-MR manganites is an exceptionally strong uniaxial 
anisotropy of the Co spins, which dramatically simplifies the possible spin arrangements. This spin 
anisotropy together with the possibility of continuous ambipolar doping turn GdBaCo2 05+a; into a 
model system for studying the competing magnetic interactions, nanoscopic phase separation and 
accompanying magnetoresistance phenomena. 

PACS numbers: 72.80.Ga, 72.20.Pa, 75.30.Kz, 75.47.De 



I. INTRODUCTION 

Since the discovery of the high-T c superconductivity 
(HTSC) in cuprates and shortly after of the colossal mag- 
netoresistance (CMR) in manganites, a great deal of ex- 
perimental and theoretical efforts have been made to clar- 
ify the nature of these phenomena. The research has soon 
revealed that the unusual behavior of cuprates and man- 
ganites is not limited to HTSC and CMR: These com- 
pounds based on seemingly simple metal-oxygen planes 
turn out to possess very rich phase diagrams, originat- 
ing from strong electron correlations and involving spin, 
charge, orbital, and lattice degrees of freedomi 2 * 3 . In 
particular, the strong electron correlations, which pre- 
vent the electrons in partially filled bands from form- 
ing conventional itinerant Bloch states, make these sys- 
tems prone to nanoscopic phase separation and self- 
organization of electrons into various superstructures. 
The role of this electron self-organization still remains 
controversial. It is often argued, for example, that the 
HTSC and CMR would never be possible in a homoge- 
neous system, and it is the nanoscopic mixture of phases 
that stays behind these novel phenomenal^ 

Ironically, the complexity of manganese and copper ox- 
ides that brings about all the fascinating physics also 
makes these compounds very difficult for understand- 
ing. The research thus naturally expanded towards 
other transition-metal oxides, both because the compar- 



ison of different systems could give a clue to the be- 
havior of cuprates and manganites, and because those 
less studied systems might be interesting in their own 
right. Such exploration has indeed been proven to 
be very fruitful, resulting, for instance, in the discov- 
ery of unconventional superconductivity in a layered- 
perovskite ruthenium oxide^ and recently in a layered 
cobalt oxide £ A study of nickel and cobalt oxides has 
also revealed the spin/charge ordering phenomena 6 - 7 - 8 9 
and nanoscopic phase separationiSiii closely resem- 
bling those in cuprates and manganites; these obser- 
vations confirm that the charge ordering is indeed a 
generic feature of strongly correlated electrons. By 
now, probably the most rich and intriguing behavior 
has been found in cobalt oxides, which ranges from gi- 
ant MR 12 i 13 i 14 i 15 i 16 i 17 and large thermoelectric power at- 
tributed to strong electron correlationsi&iSiSSiSi to un- 
conventional superconductivity. 5 Apparently, the cobalt 
oxides, which are still much less studied than cuprates 
or manganites, are meant to become the next primary 
field in investigations of the strongly correlated electron 
systems. 

Unlike cuprates and manganites, the layered cobalt 
oxides have two substantially different crystallographic 
types: the layered perovskites derived from sgwre-latticc 
C0O2 planes, similar to HTSC and CMR compounds, 
and compounds like Naa;Co02 derived from triangular- 
lattice C0O2 planes. In both cases, the C0O2 planes can 
be doped with charge carriers over a remarkably wide 



range so that the effective valence of Co ions varies from 
Co 2+ to Co 4+ . 1 In other words, the doping level ranges 
from one electron to one hole per Co ion, if the Co 3+ 
state with even number of electrons is taken as the "par- 
ent" state. Empirically, the square-lattice and triangular- 
lattice systems behave quite differently. For example, 
the square-lattice cobalt oxides, such as La^aSr^CoO^ 
Bi 2 Sr 2 Co06+5, and RBaCo 2 05 +a; (where R is a rare 
earth element), are usually reported to be non- 
metallic, 8 ! 11 ! 12 . 13 ! 14 ! 15 ! 16 . 17 ! 22 ! 23 ! 24 ! 25 with the exception of 
heavily hole-doped cubic perovskites La 1 _ 2; Sr a; Co03i 1 ^ 
In contrast, the triangular-lattice cobaltites often ap- 
pear to be fairly good metals ^S^SL with the hydrated 
Naa;Co02 being even a superconductor^ This difference 
in behavior may come in part from a disparity in the 
doping level, given that the latter compounds were usu- 
ally studied in a more highly hole-doped region. Never- 
theless, the surprisingly robust non-metallic state in the 
square-lattice cobalt oxidesi 1 ^ 5 . clearly points to a more 
fundamental source, which is presumably a very strong 
tendency to charge ordering. In the triangular-lattice 
compounds such tendency might be considerably weaker. 
Indeed, the charge ordering usually gains support from 
the antiferromagnetic (AF) exchange, while in a triangu- 
lar lattice the AF spin interactions should inevitably be 
frustrated. Thus far, however, non of these systems has 
been systematically investigated over the entire doping 
range to reveal a coherent picture of doped cobalt-oxygen 
planes. 

In this work, we undertake a systematic study of mag- 
netic and transport properties of a square-lattice cobalt 
oxide over a wide doping range. For this purpose, we 
have selected the RBaCo 2 05+ :E compounds (where R is 
a rare earth element), which have already attracted a 
lot of attention owing to such fascinating features as 
the spin-state and metal-insulator transitions, charge 
and orbital ordering phenomena, and giant magne- 
toresistance (GMR)i 8 i 9 i 12 i 13 i 14 i 15 i 16 i 17 i 22 i 23 i 24 i 28 i 29 i 30 The 
RBaCo 2 05 +I compounds possess a layered crystal struc- 
ture which consists of square-lattice layers [Co0 2 ]-[BaO]- 
[Co0 2 ]-[RO X ] stacked consecutively along the c axis (Fig. 
1) - a so-called "112"-type structured This structure 
is derived from a simple cubic perovskite Ri^Ba^MOa 
(where M is a transition metal), but in contrast to the 
latter, the rare-earth and alkali-earth ions are located in 
their individual layers instead of being randomly mixed. 

What makes RBaCo 2 05 +a: compounds particularly at- 
tractive for our study is a large variability of the oxygen 
content: By changing the annealing conditions, one can 
modify the oxygen concentration in the rare-earth [ROJ 
planes (Fig. 1) in a wide range, < x < 1 (x depends 
also on the size of R 3+ ion)M£& In turn, the oxygen con- 
tent controls the nominal valence of Co ions, which varies 
from 2.5+ to 3.5+, passing through 3+ ("parent" state) 
at x = 0.5. Ordinarily, experimental investigations of 
the T — x phase diagrams of solids are very time consum- 
ing, because they require to grow many single crystals 
of different compositions. In contrast, the composition 
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FIG. 1: Schematic picture of the RBaCc^Os+a, structure for 
x = 0.5. It is worth noting that the oxygen ions in [RO x ] 
layers demonstrate a strong tendency to ordering; for x = 
0.5, for example, they form alternating filled and empty rows 
running along the a axis ^2*—*— 



of already grown RBaCo 2 05+a; crystals can be modified 
by annealing at various temperatures and oxygen partial 
pressures, which may allow one to span the entire phase 
diagram using one and the same single crystal. Further- 
more, by varying the oxygen content one can tune the 
doping level very smoothly, which gives a great advantage 
in studying the critical regions of the phase diagram. 

Upon choosing a compound from the RBaCo^s+z 
group, one would prefer to have a non-magnetic R ion, 
such as Y, La, or Lu, to avoid additional complica- 
tion coming from the rare-earth magnetism. Unfortu- 
nately, the growth of single crystals with these non- 
magnetic elements turns out to be virtually impossible. 31 
Therefore, we have selected the GdBaCo 2 05+a; com- 
pound: Gd 3+ , being a 4I 7 ion with zero orbital mo- 
ment, is known to show rather simple magnetic be- 
havior in transition-metal compounds, making predom- 
inantly paramagnetic contribution, which can be easily 
subtracted from the overall magnetization. Also, ow- 
ing to the intermediate size of Gd ion in the series of 
rare-earth elements, GdBaCo 2 05+a; allows a fairly wide 
range of available oxygen concentration. Recently, we 
succeeded in growing high-quality GdBaCo 2 05 +:E single 
crystals using the floating-zone (FZ) technique, and stud- 
ied the magnetic and transport properties of the parent, 
x = 0.50, composition^ 

Here, we present a systematic data on the evolution 
of transport, magnetic, and thermoelectric properties of 
well-characterized GdBaCo^s+a, single crystals over the 
entire doping range available for this compound, namely, 
< x < 0.77. This layered cobalt oxide turns out to 
be a really fascinating filling-control system which allows 
a continuous ambipolar doping: We have developed a 
technique that provides an easy and precise tuning of 
the oxygen content in GdBaCo 2 05+a; single crystals and 
ceramics, and succeeded in doping the parent semicon- 
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ductor (x — 0.50) with electrons (x < 0.50) or holes 
(x > 0.50) with steps that could be as small as 0.001 
per Co ion (Ax ~ 0.001). As a result, we could ob- 
serve spectacular singularities in the transport proper- 
ties upon approaching the undoped state, could compare 
the motion of doped holes and electrons, and study the 
impact of a small amount of doped carriers on the compe- 
tition between the ferromagnetic and antiferromagnetic 
spin ordering in this compound. Our study has revealed 
a rich phase diagram for this layered cobalt oxide, which 
is found to include regions of an intriguing nanoscopic 
phase separation over virtually the entire doping range. 

This paper is organized as follows. In Section II, we 
describe the growth of high-quality GdBaCo 2 05 +;E sin- 
gle crystals by the floating-zone technique, the method 
of modifying their oxygen content, and the detwinning 
technique used to obtain single-domain orthorhombic 
crystals. The details of magnetic and transport mea- 
surements are also presented in Section II. Section III 
starts with a brief summary of the crystal structure 
of GdBaCo 2 5 + x over the oxygen-concentration range 
< x < 0.77, which is followed with the experimental 
results on transport, thermoelectric, and magnetic prop- 
erties of GdBaCo 2 5+2: , studied in this wide range of 
oxygen concentrations. The obtained data are used to 
establish an empirical phase diagram presented at the 
end of this section. The implications of our observations 
are discussed in Section IV. Based on the obtained ex- 
perimental results, first we propose a magnetic and elec- 
tronic structure of the parent GdBaCo205.so compound, 
giving our explanation of its transport behavior and the 
origin of GMR; then we discuss the effects of doping in 
GdBaCo 2 05 +a; ; and, finally, suggest an overall electronic 
phase diagram for this compound. Section V summarizes 
our findings. 



II. EXPERIMENTAL DETAILS 
A. Growth of GdBaC^Os+x crystals 

We have grown high-quality GdBaCo 2 05 +a; single 
crystals by the floating-zone technique, using an infrared 
image furnace with two halogen lamps and double ellip- 
soidal mirrors (NEC Machinery SC E-15HD). A poly- 
crystalline feed rod for the crystal growth was prepared 
by the solid-state reaction of Gd 2 03, BaC03, and CoO 
dried powders: The mixture was successively calcined 
at 850°C, 900°C, 950°C, and 1000°C, each time for 20 
hours, with careful regrinding after each sintering. Then 
the obtained homogeneous single-phase GdBaCo 2 Os +x 
powder was isostatically pressed (at ~ 70 MPa) to form 
a rod with typical dimensions of 7 mm<f> x 100 mm. Fi- 
nally, the feed rod was annealed at 1200°C in air to make 
it dense and hard. 

The crystal growth was performed in a flow of dried 
air, at a constant rate of 0.5 mm/h. We found that 
such rather slow growth rate was essential for obtain- 
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FIG. 2: (Left) X-ray Bragg's peaks for GdBaO^Os+z (x = 
0.5) crystals, which demonstrate the unit-cell doubling along 
the 6 and the c axis (each peak has CuKai and CuKq2 con- 
tributions to the diffraction pattern). For convenience, the 
peak intensity is multiplied by a factor indicated near each 
peak. (Right) X-ray (040) rocking curve. 



ing the ordered "112" crystal structure, and for reduc- 
ing the number of domains in the resulting crystalline 
rod; higher rates inevitably caused multi-domain crystal 
growth. Using optical microscopy and Laue x-ray back- 
reflection control, we selected single-domain parts of the 
grown crystal rod and cut them into parallelepiped sam- 
ples suitable for structural, transport and magnetization 
measurements. All the samples' faces were carefully pol- 
ished and adjusted to the crystallographic planes with a 
1° accuracy. 

The homogeneity and stoichiometry of the obtained 
GdBaCo 2 05+ rr crystals were analyzed by the electron- 
probe microanalysis (EPMA) and inductively-coupled 
plasma (ICP) spectroscopy which confirmed that the ac- 
tual cation composition was uniform, corresponding to 
the nominal 1:1:2 ratio within the experimental accuracy. 
Another important issue is whether Gd and Ba are well 
ordered in the lattice; in fact, we found that large rare- 
earth ions, such as La and Pr, easily mixed with Ba, 
resulting in a disordered cubic phase Ro.sBao.sCoOa-^. 
In the case of GdBaCo 2 05 +a; crystals, however, the x-ray 
diffraction data demonstrate that Gd and Ba are indeed 
well ordered into consecutive (001) layers, which results 
in the doubling of the unit cell along the c axis (Fig. 2). 
Moreover, for the oxygen concentration x « 0.5, the oxy- 
gen ions are found to form alternating filled and empty 
chains running along the a axis; this brings about the 
doubling of the unit cell along the b axis (Fig. 2). Appar- 
ently, such long-range oxygen ordering would hardly be 
possible if considerable amount of Ba were substituting 
Gd in GdOa, layers. Figure 2 shows also a typical x-ray 
rocking curve (040) for GdBaCo 2 05.5, which has a full- 
width-at-half-maximum (FWHM) of less than 0.1°, indi- 
cating that our crystals have few macroscopic defects. An 
additional evidence for the macroscopic crystallographic 
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perfection is the ease with which the crystals could be 
cleaved, especially along the {001} planes, exposing per- 
fectly flat, shiny surfaces. 



B. Tuning the oxygen content in GdBaC^Os+z 

1. Equilibrium oxygen concentration as a function of 
temperature and oxygen partial pressure 

Owing to the layered "112" crystal structure, the oxy- 
gen content in RBaC^Os+z can be varied within a wide 
range < x < 1, where the oxygen vacancies are lo- 
cated predominantly in the rare-earth RO^ planesJ^ The 
most convenient way to modify the oxygen stoichiometry 
is the high-temperature annealing, yet it requires a de- 
tailed knowledge of how the equilibrium oxygen content 
depends on the temperature and oxygen partial pressure. 
Besides, one needs to know the kinetics of the oxygen ex- 
change: The annealing temperature should be carefully 
chosen so that the oxygen exchange is quick enough for 
the equilibrium state to be reached in a reasonable time, 
but still slow enough to avoid unwanted oxygen absorp- 
tion during the subsequent cooling of the sample. 

In order to establish such dependences, we have per- 
formed special sets of annealings, systematically varying 
temperature, oxygen partial pressure P02, and anneal- 
ing time, as exemplified by the data shown in Fig. 3. 
By measuring the weight change of ~ 1-g polycrystalline 
samples, as well as large single crystals, with a 0.1-/ig res- 
olution, we could evaluate the change in the oxygen con- 
tent Ax with an accuracy better than 0.001. Although 
all the relative oxygen variations are measured very accu- 
rately, there would still be a large uncertainty in the abso- 
lute values of x, unless we pin this relative scale to the ab- 
solute oxygen content at least at one point. Fortunately, 
in GdBaCo 2 05 +2; there are two peculiar oxygen concen- 
trations that allow the absolute x scale to be establish 
unambiguously. First, in GdBaC^Os+a; crystals, similar 
to other RBaC^Os+a; compounds with small rare-earth 
ionS ) 8 i 9 i 28 the oxygen content can be reduced down to 
x = by annealing in vacuum, or inert atmosphere. One 
can naturally expect the x = composition to be stable 
over rather broad range of parameters, since in this phase 
all the weakly-bound oxygen is removed from CdO^ lay- 
ers, while the strongly-bound oxygen in BaO and C0O2 
layers is still intact. Indeed, we observed that x satu- 
rates, approaching some lowest level [see Fig. 3(b)] as 
the oxygen partial pressure is reduced and the annealing 
temperature is increased, and we attributed this satu- 
ration value to x = 0. Much more sensitive calibration 
can be done at the x = 0.5 point, which turns out to be 
critical for GdBaCo 2 05 +2; : Upon crossing this point, the 
mixed valence composition of Co 2+ /Co 3+ ions turns into 
the Co 3+ /Co 4+ one, and consequently the type of charge 
carriers switches abruptly. As we will show below, sev- 
eral physical properties of GdBaCo205+a; exhibit a very 
sharp singularity in their x dependences near x = 0.5, 
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FIG. 3: The equilibrium oxygen concentration a; oq in 
GdBaCo2 05+!. (a) Temperature dependences of a; cq for sev- 
eral values of the oxygen partial pressure, (b) Dependences 
of x eq on the oxygen partial pressure at several temperatures. 

changing remarkably with a minute (by merely 0.001) 
modification of x. The calibrations that used the x = 
and x — 0.5 points gave the same results, which allows 
us to be confident in absolute values of x in our samples. 

Temperature dependences of the equilibrium oxygen 
content x in GdBaCo 2 05 +2; for different oxygen partial 
pressures are shown in Fig. 3(a), which demonstrates 
that a large variation of x can be achieved by relatively 
simple means: by annealing in a flow of oxygen and ar- 
gon mixed in different proportions. Note that the shown 
data correspond to the equilibrium state, and that all the 
changes in oxygen concentration are therefore completely 
reversible. Our measurements have shown that the equi- 
librium x value at a given temperature is roughly propor- 
tional to the logarithm of the oxygen partial pressure, as 
long as x is not too close to zero [Fig. 3(b)]. Therefore, 
in order to obtain the lowest oxygen concentration x = 
we annealed samples in a flow of argon or helium, and 
used strongly diluted (down to 10 ppm) oxygen- argon 
mixtures to access the low- a; range, although precise tun- 
ing of x in this case becomes technically more difficult. 
In the opposite limit, the samples with oxygen concen- 
trations up to x = 0.77 were prepared by annealing in 
oxygen at high pressures (up to 70 MPa). 



2. Oxygen intercalation kinetics 

As we have mentioned above, there is one more point 
to be concerned about, upon choosing the proper an- 
nealing conditions, besides the equilibrium x value: it 
is the oxygen-intercalation kinetics. If the oxygen diffu- 
sion is too slow, one will inevitably end up with a crystal 
having a large composition gradient; if, in contrast, the 
oxygen uptake is too fast, it will be difficult to preserve 
the achieved state even by very fast cooling. By measur- 
ing how the sample's mass evolves with time t, we have 
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FIG. 4: Kinetics of the oxygen intercalation into GdBaCo2 05+ I . (a) Comparison of the oxygen uptake in a ceramic sample 
(with an average grain size of ~ 10 (im) and a single crystal (with sizes in the ab plane of 500 fim x 400 fim) measured in 
the oxygen flow at 250° C. Solid lines are the result of simulation, see text, (b) The oxygen uptake in a single crystal with 
sizes in the ab plane of 1000 fj,m x 600 /im measured in the oxygen flow at different temperatures (solid lines are the result of 
simulation), (c) Dependence of the oxygen-uptake kinetics in a ceramic samples at 600° C on the oxygen partial pressure. 



determined the time dependences of the average oxygen 
content x{t) for different annealing conditions, as shown 
in Fig. 4. Apparently, the proper duration of the heat 
treatment necessary to obtain a homogeneous oxygen dis- 
tribution should be several times longer than the charac- 
teristic time r of the oxygen exchange process at a given 
temperature and an oxygen partial pressure. We have 
found that for ceramic samples the kinetics of oxygen in- 
tercalation in the entire studied ranges of temperature 
and P02 follows a simple exponential law, 

x(t) = - [xoo - a; ]e _t/T , (1) 

where xq and the initial and the equilibrium oxy- 

gen contents, r is a time constant that depends on tem- 
perature and Po2- Such simple exponential dependence 
can be naturally attributed to the oxygen exchange lim- 
ited by a surface energy barrier. The time constant r 
turns out to be less than 3 hours at temperatures down 
to 250° C (at 1 bar oxygen pressure), so that the equilib- 
rium state in ceramic samples can be easily achieved by 
one-day annealing [Fig. 4(a)]. 

In the case of single crystals, however, the oxygen inter- 
calation deviates from the simple exponential behavior, 
and goes at a noticeably slower rate, as illustrated in Fig. 
4(a). We have found that for ~ 500 fim sized crystals, 
two days of annealing at 250° C is far from being enough 
to reach the equilibrium oxygen concentration [Fig. 4(a)]; 
such annealing turns out to be insufficient even for ~ 100 
/jm samples. This difference in the kinetics of oxygen ab- 
sorption is apparently caused by an additional limitation 
imposed by the bulk oxygen diffusion: While the crystal 
surface layer is readily filled with oxygen, it takes much 
more time for oxygen to diffuse towards the inner part 
of the crystal. For rectangular-shape crystals the oxygen 
exchange kinetics has an analytical solutionjiSS which 



depends on the surface exchange coefficient K, character- 
izing the oxygen exchange at the interface between the 
gas and the solid, and the chemical diffusion coefficient 
D, besides the sample size. These parameters can be ex- 
tracted by fitting theoretical curves to experimental data, 
as exemplified in Figs. 4(a) and 4(b). For example, the 
data shown in Fig. 4(a) provide the diffusion-coefficient 
value of D = 3 x 10~ 8 cm 2 /s at 250 o C^4 With increasing 
temperature, the diffusion coefficient grows rapidly and 
so does the rate of the oxygen uptake, making it possi- 
ble to achieve a homogeneous oxygen distribution in ~ 1 
mm sized crystal already after several hours annealing at 
400°C [Fig. 4(b)]. Note that the oxygen kinetics does not 
show any detectable difference upon varying the sample 
size along the c axis, implying an essentially 2D character 
of the oxygen diffusion typical for layered oxides^ 

We should note that the diffusion coefficient found in 
GdBaCo205+ K is remarkable in its own right, being un- 
usually large for such low temperatures; the oxygen dif- 
fusivity appears to be comparable with that in best su- 
perionic conductors^ One additional implication of the 
high oxygen mobility, which is important to the present 
study, is that the oxygen ions may be capable of rearrang- 
ing even at room temperature or below, and one should 
expect the oxygen to form ordered superstructures or to 
participate in mesoscopic phase separation. 

After the homogeneous oxygen distribution in a crys- 
tal is reached, an important issue to be concern about 
is the cooling procedure, which must be as fast as pos- 
sible to avoid any further oxygen uptake. To obtain the 
most homogeneous samples we used the following rules 
upon selecting the annealing conditions. The anneal- 
ing temperature and the partial oxygen pressure were 
selected so that: i) they provided the required oxygen 
content in the crystal, and ii) the annealing time nec- 
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essary to reach the equilibrium oxygen distribution was 
in the range from several hours to several days. After 
the annealing was completed, the samples were quickly 
quenched to room temperature without changing the 
atmosphere. Upon quenching, the crystal temperature 
dropped by 100 — 200°C in several seconds, guaranteeing 
the blocking of any further oxygen exchange. The huge 
difference in the time scales - hours and days for modify- 
ing the oxygen concentration and seconds for quenching 
- ensured us that the unwanted oxygen uptake could af- 
fect no more than a fraction of percent of the crystal's 
volume. An important point to be specially emphasized 
is that the oxygen exchange at the interface between the 
gas and the solid is strongly suppressed at low oxygen 
partial pressures. As illustrated in Fig. 4(c), by reducing 
the oxygen partial pressure one can arrange a slow oxygen 
uptake in a ceramic sample even at 600° C. This feature 
turns out to be very useful in obtaining homogeneous 
samples within the low oxygen concentration range. 



C. Detwinning of crystals 

The oxygen ordering is indeed observed in 
GdBaCo 2 05 +a; at x « 0.5, where the oxygen ions 
in GdOz planes order into alternating filled and empty 
chains, causing a tetragonal-to-orthorhombic (T-O) 
transition and doubling of the unit cell along the b axis 
(Fig. 2). Usually this T-0 transition is accompanied 
by heavy twinning of crystals that mixes the a and b 
orthorhombic axes; one therefore needs to perform a de- 
twinning procedure to get a single-domain orthorhombic 
crystal for studying the in-plane anisotropy. To detwin 
crystals, we slowly cooled them under a uniaxial pres- 
sure of ~ 0.15 GPa from 260°C, using a polarized-light 
optical microscope to control the twin removing^ We 
should note, however, that the GdBaCo 2 05 +a; crystals 
are very fragile, and it was virtually impossible to 
complete the detwinning procedure without having them 
cleaved into two or more pieces. These pieces were fairly 
well detwinned (according to x-ray measurements, the 
remaining fraction of misoriented domains usually did 
not exceed 4-5%) and suitable for magnetization mea- 
surements. Thus far, however, we have not succeeded in 
detwinning the crystals with already prepared electrical 
contacts (see below), which are necessary for transport 
measurements. 



D. Details of measurements 

To characterize the physical properties of 
GdBaCo205+a;, we have carried out resistivity, mag- 
netoresistance (MR), Hall, thermopower, and magneti- 
zation measurements within the 2-400 K temperature 
range. 

The in-plane (p a b) and out-of-plane (p c ) resistivity was 
measured using a standard ac four-probe method. Good 



electric contacts were obtained by drawing with a gold 
paint on polished crystal surfaces, and subsequent heat 
treatment. For current contacts, the whole area of two 
opposing side faces was painted with gold to ensure a 
uniform current flow through the sample. In turn, the 
voltage contacts were made narrow (~ 50 /im) in order to 
minimize the uncertainty in absolute values of the resis- 
tivity. It is important to note that the gold contacts were 
prepared before all the heat treatments that were used 
to vary the oxygen content. After the required oxygen 
concentration was set by annealing, thin gold wires were 
attached to the contact pads using a room-temperature- 
drying silver paste (DuPont 4922), which electrically and 
mechanically bound the wire to the sample. As the mea- 
surements were completed, the wires were removed and 
the crystal was reannealed to get the next x value. 

The MR measurements were done either by sweeping 
the magnetic field between ±14 T at fixed temperatures 
stabilized by a capacitance sensor with an accuracy of 
~1 mK, 37 or by sweeping the temperature at a fixed 
magnetic field. Both Ap c / p c and Ap a b/p a i, have been 
measured for H || ab and H || c. 

The Hall resistivity was measured using a standard 
six-probe technique by sweeping the magnetic field H || c 
to both plus and minus polarities at fixed temperatures; 
the electric current was always along the a6-plane. After 
the symmetric (MR) contribution coming from a slight 
misplacement of the contacts was subtracted from the 
raw data, the Hall resistivity appeared to be perfectly 
linear in magnetic field for all measured temperatures, 
implying that the anomalous Hall effect was negligible. 

In order to determine the thermoelectric power (See- 
beck coefficient) S, we generated a slowly oscillating ther- 
mal gradient along the sample (within ~ 1 K), and mea- 
sured an induced periodic voltage. This allowed us to 
get rid of thermoelectric contributions generated in the 
remaining circuit. A chromel-constantan thermocouple 
employed for measuring the thermal gradient was at- 
tached to the heat source and to the heat sink, which 
were fixed to the sample by a silver paste. To obtain 
the absolute value of the thermoelectric power, a contri- 
bution from the gold wires (~ 2 pV/K) used as output 
leads was subtracted. Thermopower measurements were 
mostly performed on long (> 2 mm) samples by applying 
a temperature gradient along the a or the b axis. 

Magnetization measurements were carried out using a 
SQUID magnetometer at fields up to 7 T applied along 
one of the crystallographic axis. Measurement modes 
included taking data upon heating the sample after it 
was cooled down to 2 K in zero field (ZFC), upon cooling 
the sample in magnetic field from 400 K (FC), and upon 
sweeping the field at fixed temperatures. Throughout 
this paper, the magnetization coming from Co ions is 
determined by subtracting the contribution of Gd ions, 
assuming their ideal paramagnetic (PM) behavior with 
total spin S = 7/2; the latter is a good approximation 
since no ordering of Gd 3+ moments in any of the samples 
is detected down to 1.7 K. 
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III. RESULTS 

A. Crystal structure 

Figure 5 shows the evolution of the GdBaCo205-|_ a ; 
crystal structure as a function of the oxygen concentra- 
tion x. We can clearly distinguish three different compo- 
sition regions: 

• < x < 0.45, the system keeps a macroscopically 
tetragonal structure, where the unit cell smoothly 
expands in the c direction and shrinks in the in- 
plane directions with increasing x; 

• 0.45 < x < 0.60, the oxygen ions order into alter- 
nating filled and empty chains running along the 
a axis, which results in the orthorhombic structure 
and in the doubling of the unit cell along the b axis 
(see Fig. 2); 

• x > 0.60, the system evolves towards a macroscop- 
ically tetragonal symmetry, though in quite a com- 
plicated way. For crystals located at the lower edge 
of this composition range, we observed at the same 
time the x-ray diffraction peaks corresponding to 
the orthorhombic structure with the unit cell dou- 
bling, and those related to the tetragonal structure. 
This two-phase state is clearly intrinsic and unre- 
lated to a macroscopic oxygen-concentration gradi- 
ent that would emerge if crystals were improperly 
annealed. Indeed, the signs of the two-phase state 
were reproducibly observed only in a quite narrow 
range of x and disappeared as the oxygen content 
was further increased. At the highest achieved oxy- 
gen concentration, we could distinguish only the 
tetragonal, albeit somewhat broadened, diffraction 
peaks. 
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FIG. 5: Room-temperature parameters of the GdBaC^Os+a, 
unit cell as a function of the oxygen content. 



We should keep in mind, however, that this dividing is 
based only on the macroscopic symmetry, and thus it 
should not be taken too literally; the behavior of the 
local structure may in fact be much more tricky. For 
example, it is quite likely that oxygen ordering or meso- 
scopic phase separation takes place at small x values as 
well, but the ordered domains are too small to be dis- 
cernible by conventional x-ray diffraction. The same is 
also true for high oxygen contents where one may expect 
fairly ordered phases at x « 2/3; 3/4, as well as various 
mesoscopic phase mixtures. 



B. Resistivity 

In perovskite compounds, the electron-band filling can 
be modified by a partial substitution of cations with el- 
ements having a different valence, as in the case of the 
Lai-^Sr^Af O3 compounds (M is a transition metal), and 
by introducing cation or anion vacancies and interstitials 
as in V2- y O-3, LaTi03 + a, and YBa2Cu3G 1 6+x'^ At in- 
teger filling (integer number of electrons per unit cell) 
these compounds are usually band or Mott insulators, 
yet a metallic state often emerges upon changing the fill- 
ing level, that is, when electron or holes are doped into 
the parent insulator. 

In GdBaCo205+a;, it is obviously the x = 0.5 com- 
position that is the parent compound, where all the Co 
ions are nominally in the 3+ state. The limiting phases 
x = and x = 1 should correspond to the 1:1 mixtures 
Co 2+ /Co 3+ and Co 3+ /Co 4+ , respectively, or in other 
words to the doping levels of 0.5 electrons and 0.5 holes 
per Co ion. It would not be surprising, therefore, if the 
x = 0.5 composition were insulating, and a metallic be- 
havior were emerging when the oxygen content deviates 
from x = 0.5 towards lower or higher values. The ac- 
tual behavior of GdBaCo205+ a: , however, turns out to 
be more complicated: In contrast to the naive expec- 
tations, it never becomes a true metal (Fig. 6), even 
though we change the doping level in a very broad range 
from 0.5 electrons per Co (x = 0) up to 0.27 holes per 
Co (x — 0.77)m& At low temperatures, the in-plane re- 
sistivity Pab exhibits an insulating behavior throughout 
the whole accessible range of x. Moreover, by looking at 
the low-temperature region (T < 150 K) in Fig. 6, one 
may easily find out that the resistivity even increases as 
the oxygen content is reduced below 0.5 ("electron dop- 
ing"); clearly, the observed evolution differs remarkably 
from what one usually expects to see upon doping an 
insulator. 

Although GdBaCo205+a; never behaves as a normal 
metal, the temperature dependence of its resistivity does 
not follow that of a simple insulator or semiconductor 
either. As the temperature increases above 200 — 250 K, 
the in-plane resistivity p a b(T) shows a gradual crossover 
(for x < 0.44) or a sharp transition at T w 360 K 
(for x close to 0.5) into a metal-like state, see inset of 
Fig. 6. This sharp "metal-insulator" transition (MIT) 
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FIG. 6: In-plane resistivity p a b{T) of GdBaCo^Os+rc crystals 
with different oxygen concentrations. Inset: expanded view 
of p a b{T) in the vicinity of the "metal-insulator" transition. 
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eral temperatures. 



at x » 0.5 is similar to what has been reported for air- 
sintered RBaCo 2 05 +a; ceramicsi 12 i 13 i 14 i 15 Although the 
resistivity p ab ~ 400 — 600 pfl cm on the metal side of 
the MIT is still relatively large, in transition-metal ox- 
ides such resistivity values are actually more consistent 
with a metallic, rather than a hopping transport^ As 
can be seen in the inset of Fig. 6, the metal-insulator 
transition is the sharpest and takes place at the highest 
temperature Tmit ~ 364 K for x = 0.50, while any de- 
viation from this stoichiometry blurs the MIT and shifts 
it to lower temperatures. This is particularly clear for 
high oxygen contents: Tmit is reduced to ~ 325 K for 
x = 0.65, and no distinct transition can be found for 



x > 0.70. In addition to the MIT, a clear, albeit small, 
kink is seen on the p(T) curves in the temperature range 
100-260 K, which is related to a magnetic transition as 
will be discussed below. 

To make it more transparent how the resistivity in 
GdBaCo 2 05 +x evolves with oxygen content, in Fig. 7 
we plot pab at several temperatures as a function of x. 
Regardless of the temperature, the ^-dependence of re- 
sistivity splits into roughly the same three composition 
regions as were observed in the behavior of the crystal 
structure (Fig. 5). For both low (0 < x < 0.45) and 
high (x > 0.7) oxygen concentrations, the resistivity ap- 
pears to be almost rr-independent, while in the region 
0.45 < x < 0.7 it changes rather steeply, by orders of 
magnitude. This step-like resistivity drop, taking place 
upon going from "electron-doped" to "hole-doped" com- 
positions, tends to hide the singularity expected for the 
parent compound. Only very detailed data collected with 
a step of Sx < 0.01 make this singularity discernible on 
the low-temperature p a t—x curves, where it is manifested 
as a narrow peak in the vicinity of x = 0.50. 

In fact, the most unusual and intriguing feature in Fig. 
7 is the asymmetry with respect to the oxygen concentra- 
tion x = 0.50. The resistivity evolution at x > 0.50 looks 
rather conventional: Indeed, doping of ~ 2.5% of holes 
(x = 0.525), which turns Co-ions into a mixed 3+/4+ va- 
lence state, dramatically improves the conductivity. As 
can be seen in Fig. 7, the in-plane resistivity at 100 
K drops by more than two orders of magnitude as x 
increases from « 0.50 to 0.525, and by three orders of 
magnitude more upon further increasing x. In contrast, 
upon electron doping (reducing x below 0.50), an initial 
decrease of the resistivity almost immediately turns into 
a resistivity growth. Moreover, in a wide region x < 0.44 
both the absolute value of the resistivity and its tempera- 
ture dependence are virtually independent of the oxygen 
content (Figs. 6, 7). Such insensitivity of the conduc- 
tivity to doping is very unusual, and may only be pos- 
sible if the electrons released upon removing the oxygen 
do not participate in the charge transport. This clearly 
speaks against the "rigid-band" picture, where electrons 
are expected to fill the lowest-energy unoccupied states. 
Apparently, the electron doping is accompanied by de- 
veloping of microscopic insulating states or mesoscopic 
regions, where the electrons are immediately trapped. 
At this point it is useful to remind that the insulating 
nature of the limiting x = composition has been al- 
ready understood as coming from the charge ordering of 
electrons in C0O2 planes into a unidirectional charge- 
density wavers It is therefore reasonable to suggest that 
the robust insulating behavior observed in the wide range 
< x < 0.44 is also associated with some kind of charge 
ordering among the doped electrons. 

It is interesting to analyze the low-temperature behav- 
ior of the resistivity, which can give information on the 
transport mechanisms operating in the system. Data fit- 
ting has shown that a crystal with the precisely tuned 
x = 0.50 composition exhibits a simple activation behav- 
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x = 0.45 - 0.525; and (c) Efros-Shklovskii behavior for x = 0.165 - 0.440. 



ior p ab oc exp(A/T) with A « 70 meV [Fig. 8(a)], which 
is well consistent with the understanding of this com- 
position as a parent insulator (narrow-gap semiconduc- 
tor). However, when the oxygen content deviates some- 
what from x = 0.50, the resistivity behavior almost im- 
mediately switches into the 3D variable-range hopping 
(VRH) modej^a Pab oc exp[(T /r) 1 / 4 ], as shown in Fig. 
8(b). This type of conduction is typical for disordered 
systems where the charge carriers move by hopping be- 
tween localized electronic states. The formation of such 
localized states may be rather simply conceived of in the 
following way: Initially, at the oxygen composition of 
x = 0.50, GdBaCo 2 05 +2; possesses a well-ordered crystal 
structure, where the oxygen ions form perfect filled and 
empty chains alternating along the b axis (Fig. 1). When 
the oxygen concentration deviates from x = 0.50, this re- 
sults either in vacancies emerging in the filled chains, or 
in oxygen ions that go into the empty chains. While these 
oxygen defects inevitably generate electrons or holes in 
the C0O2 planes, they also produce a poorly screened 
Coulomb potential that may well localize the generated 
carriers, so that some of the adjacent Co ions acquire 
the Co 2+ or Co 4+ state. Then the conductivity oc- 
curs through hopping (tunneling) motion of such local- 
ized Co 2+ or Co 4+ states. As can be seen in Fig. 8(b), 
the slope of the resistivity curves monotonically decreases 
with increasing x, implying that the localization length 
of holes (Co 4+ states) is noticeably larger than that of 
electrons. 

Upon further decreasing the oxygen content below 
x « 0.45, the temperature dependence of the resistiv- 
ity becomes steeper than expected for the Mott's VRH 
regime, p a b oc exp[(To/T) 1 / 4 ], and we find that it is 



better fitted by the Efros-Shklovskii expression for the 
hopping conductivity^ p a b oc exp^Tp/T) 1 / 2 ], see Fig. 
8(c). The latter behavior is usually observed when the 
Coulomb interaction starts to play a key role in carriers 
hopping, bringing about a strong depletion in the den- 
sity of states (Coulomb gap) near the Fermi energy. The 
data in Fig. 8 imply that the Coulomb-repulsion effects 
in GdBaCo 2 05 +2; gain strength as the oxygen content is 
reduced below « 0.45. 

The unusual resistivity behavior at high oxygen con- 
centrations (Fig. 6) is also worth mentioning. For both 
the x = 0.65 and x = 0.70 crystals, p a f, smoothly in- 
creases by 2 — 3 orders of magnitude upon decreasing 
temperature, and reaches apparently non-metallic values 
of the order of 1 ficm. Nevertheless, it becomes also 
clear from Fig. 6 that, at least in the temperature range 
studied, the resistivity tends to saturate at a finite value, 
instead of diverging at T = 0; such behavior would in- 
dicate a metallic ground state, if the resistivity values 
were not that large. A possible solution of this puzzle is 
an intrinsic mesoscopic phase separation, that makes the 
carriers to move along filamentary conducting paths^i 
In this case, the zero-temperature conductivity, being de- 
termined by the topology of the metallic phase, may be 
arbitrary small and can easily violate the Mott limit 39 
that sets the minimum metallic conductivity for homo- 
geneous systems. This picture is quite plausible, since, as 
we have mentioned above, the crystal structure at these 
high oxygen concentrations actually bears signs of phase 
separation. 

Thus far we discussed the resistivity behavior along 
the C0O2 planes, which is however not the full story, 
since one might anticipate a strong electron-transport 
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FIG. 9: The resistivity anisotropy in GdBaCc^Os+a:- (a) Resistivity of GdBaCc^Os.so measured along the c axis and along the 
ab plane, (b) The resistivity anisotropy ratio p c / p a b for x = 0.165 — 0.440 and 0.500. 



anisotropy to be brought about by the layered crys- 
tal structure of GdBaCo205+ :r . Indeed, many lay- 
ered transition-metal oxides, such as high- J), cuprates, 42 
manganites^ or even cobaltites built from triangular- 
lattice C0O2 planesji2*2& exhibit huge anisotropy val- 
ues and contrasting temperature dependences of the 
in-plane and out-of-plane resistivity. In the case of 
GdBaCo205+a;, where the CdO^ layers with variable 
oxygen content are located in between the C0O2 planes 
(Fig. 1), it would be natural to expect a kind of 3D-to-2D 
transition to occur upon reducing x from 1 to 0, that is, 
as the oxygen ions binding the C0O2 planes are removed. 
However, a comparison of the in-plane and out-of-plane 
resistivity in Fig. 9 shows that this is not really the case. 
The anisotropy p c / pab indeed increases slightly as the 
oxygen concentration is reduced from x ~ 0.44 towards 
zero [Fig. 9(b)], but remains rather moderate. In fact, a 
considerable anisotropy is observed only at intermediate 
temperatures, while at both high and low temperatures 
GdBaCo205+2; tends to become virtually isotropic. It 
turns out therefore that the oxygen-depleted GdO^ layers 
do not constitute a serious obstacle for the c-axis electron 
motion. 

Somewhat different behavior of the parent 
GdBaCo2C>5.5o compound provides a clue to un- 
derstand the mechanism responsible for resistivity 
anisotropy. As can be seen in Figs. 9(a), 9(b), the 
in-plane and out-of-plane resistivities stay virtually 
indistinguishable from each other down to T « 260K, 
where the FM-AF transition takes placed Upon further 
decreasing the temperature, the curves sharply diverge 
and p c / pab grows from w 1 at T > 260K up to » 22 
at 100 K. This indicates that the charge motion in 
GdBaCo205+a; is very sensitive to the arrangement 
of Co spins, and it is the peculiar spin ordering in 



this system that is responsible for the appearance of 
conductivity anisotropy. 

C. Thermoelectric power 

Among the features that currently attract a lot of 
attention to transition-metal oxides are their unusual 
and potentially useful thermoelectric properties ^iM^^i 
The peculiar thermoelectric behavior in TM oxides is of- 
ten attributed to strong electron correlationSf22i2I though 
the picture still remains far from being clear. In this 
respect, the RBaCo205 +:r compounds, being capable of 
smoothly changing the doping level in a very broad range, 
appear to provide a good testing ground for studying the 
problem. 

The thermoelectric power S(T) measured on both 
single-crystal and ceramic GdBaCo205+a; samples^ 
throughout the available oxygen-concentration range is 
presented in Fig. 10. Following the usual approach in 
analyzing the behavior of semiconductors, we have plot- 
ted the data in the inverse-temperature scale, since in 
semiconductors the thermoelectric power is expected to 
be linear in 1 /T with a slope reflecting the activation en- 
ergy for charge carriers A s , S(T) « A±(kB/e)(A s /kBT). 
It becomes immediately clear from Fig. 10 that the only 
composition that exhibits a conventional semiconducting 
behavior is x = 0.500, with the activation energy A s w 70 
meV being in good agreement with that deduced from 
the resistivity data in Sec. IIII Bl Just a subtle deviation 
of the oxygen concentration from x = 0.500 (by merely 
0.001) qualitatively changes the thermoelectric behavior: 
S becomes almost temperature independent at w 100 K, 
indicating that the electron transport is no longer gov- 
erned by the band gap. It is worth noting that the resis- 
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FIG. 10: Thermoelectric power of GdBaCo205+a; for different 
oxygen concentrations plotted as a function of inverse temper- 
ature. 



tivity of GdBaCo205+j; also switches its behavior from 
a simple activation one into the variable-range hopping 
mode upon deviation from the parent composition (Fig. 
8). 

The salient feature in Fig. 10 is a very abrupt change 
in sign of the low-temperature thermoelectric power upon 
crossing x = 0.500; 5(100 K) jumps from w -700 /xV/K 
at x = 0.497 up to « +800 /zV/K at x = 0.501. In 
fact, this sign change unambiguously indicates that the 
type of change carriers sharply switches from electrons 
(x < 0.5) to holes (x > 0.5) without any messy inter- 
mediate state. It turns out to be possible, therefore, 
to drive the doping level in GdBaC^Os+z continuously 
across the parent insulating state (so that the Fermi level 
jumps across the gap), which is really unusual among 
transition- metal oxides. The feature that makes possible 
such continuous and smooth doping of GdBaCo205 +:E is 
the metal- insulator transition at Tmit: At temperatures 
where the oxygen concentration is modified, the gap in 
the electronic band structure is closed and, therefore, the 
chemical potential changes gradually with the changing 
of the carriers concentration. 

The metal-insulator transition at Tmit ~ 360 K has a 
spectacular manifestation in the thermoelectric proper- 
ties (Figs. 10 and 11), owing to the high sensitivity of 
5(T) to details of the electronic structure. As can be seen 
in Fig. 11, at high temperatures T > 370 K, 5 is very 
small (« —4 /zV/K), as is usually the case with ordinary 
metals, and virtually independent of temperature. More- 
over, it does not change with doping either, implying that 



the Fermi level is positioned in the middle of a wide con- 
duction band. Upon cooling the sample below 370 K, 
however, the thermoelectric power abruptly increases, in- 
dicating that the conduction band splits and a band gap 
opens near the Fermi level. Depending on the doping, 
the Fermi level gets trapped inside the newly created va- 
lence or conduction band, determining which type of car- 
riers - electrons or holes - will prevail. Immediately after 
the transition, 5 appears to be positive even in slightly 
electron-doped samples (x < 0.5), indicating that while 
both electrons and holes are contributing to charge trans- 
port at these temperatures, the holes have considerably 
higher mobility. At lower temperatures, carrier excita- 
tions over the gap become suppressed, letting only one 
type of carriers survive; consequently, at T < 250 — 280 K 
the sign of 5 is uniquely determined by the doping level 
(x < 0.5 or x > 0.5). 

One may wonder whether and how these two transport 
properties - resistivity and thermoelectric power - corre- 
late in GdBaCo205+2;. In conventional doped semicon- 
ductors, for instance, such correlation is rather straight- 
forward if one type of carriers dominates. Namely, the 
resistivity and thermoelectric power are both governed by 
the activation energy of carriers, p oc exp(A/fcsT) and 
5 « A± (fce/e)(A/fcBT), and thus are related roughly as 
5 C± (fcs/e) In p. This qualitative relation implies that 
upon doping a semiconductor, its thermoelectric power 
should decrease (in absolute value) by roughly 200 /iV/K 
per order-of-magnitude reduction in resistivity. In fact, 
this trend is not restricted to conventional semiconduc- 
tors, but holds also in many other non-metallic systems, 
including doped Mott insulators42^24S though the slope 
dS/d(\np) may differ from fc B /e » 86.2 pV/K. 

If this S-vs.-p correlation were also working in the case 
of GdBaCo205+a;, its thermoelectric power would show a 
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FIG. 12: Doping dependences of the thermoelectric power in 
GdBaC^Os+z at several temperatures. 



step-like decrease in magnitude upon increasing x, follow- 
ing the asymmetric resistivity curves (lnp vs. x) shown 
in Fig. 7. However, Fig. 12, which presents S at several 
temperatures as a function of the oxygen concentration, 
makes it immediately clear that these expectations fail. 
Instead, the doping dependence of thermoelectric power 
turns out to be fairly symmetric with respect to the par- 
ent, x — 0.5, composition: S(x) exhibits a spectacular 
divergence with approaching x = 0.5, where it reaches 
large absolute values and changes its sign upon crossing 
this peculiar doping level. It is worth noting that S(x) 
switches from being negative (electron-like) to positive 
(hole-like) within a very narrow but discernible region 
around x = 0.5 (its expanded view is given in Fig. 13); 
this slight thermally-induced smearing gives evidence for 
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FIG. 14: Temperature dependences of the thermoelectric 
power in GdBaCo2 05+ I for x = 0. 



a continuous evolution in the Co02-plane doping. When 
the oxygen concentration deviates from x — 0.5, S(x) 
smoothly reduces its magnitude, eventually approaching 
a small constant value (Fig. 12). The S(x) curves keep 
this singular and fairly symmetric behavior regardless of 
temperature, except for the high-temperature "metallic" 
region (T > 370 K), where the thermoelectric power is 
small and virtually doping independent (Fig. 12). In 
fact, the striking contrast in the doping dependences of 
two transport properties - symmetric S(x) and asym- 
metric p(x) - is one of the most intriguing feature in the 
thermoelectric behavior of GdBaC^O^+a;. 

The smooth evolution of S with decreasing oxygen con- 
tent below x = 0.5 shown in Fig. 12 breaks down, how- 
ever, upon approaching the limit of x = - another 
peculiar composition corresponding to the 50% electron 
doping. This composition is known for the long-range 
charge ordering that sets in below w 210 — 220 K, as was 
observed in YBaC^Os.o and HoBaCo2C>5.o»^ A similar 
charge ordering presumably takes place in GdBaC^Os.o 
as well, and it has a clear manifestation in the behavior of 
thermoelectric power, whose magnitude quickly increases 
below 230-240 K (Fig. 14). The long-range charge order- 
ing at 50%-electron doping and its impact on the charge 
transport, being interesting in their own right, need a 
detailed study, which however goes beyond the scope of 
the present paper. 



D. Magnetization 

1. General features 

Broadly speaking, it would be quite natural for 
cobalt oxides to possess magnetic properties more com- 
plex than those of isostructural compounds based on 
most other transition metals. The reason for this 
additional complexity is that cobalt ions in a given 
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FIG. 15: Raw magnetization M(T) data obtained for 
GdBaCo2 05+i single crystals with x = 0.25 (a) and x = 0.50 
(b) at the magnetic field H = 7 T applied parallel or per- 
pendicular to the c axis. Dashed lines show the Curie- Weiss 
contribution of Gd 3+ ions (/i e g = 7.94 /ib; 9 = K). 



valence state can have more than one allowed spin 
state; for example, Co 3+ ions can acquire a low-spin 
(LS: t% g , e°; S = 0), intermediate-spin (IS: t\ g , 
ej; S = 1), or high-spin (HS: 4 g , e 2 ; S = 
2) state, 8 ! 9 ! 10 ! 11 ! 12 . 13 ! 14 ! 15 ! 16 ! 17 ! 22 ! 23 ! 24 ! 28 ! 29 ! 30 These states 
appear to be located quite close to each other in en- 
ergy, bringing about a possibility of the spin-state tran- 
sitions/crossovers upon changing temperature or lattice 
deformation. 

Keeping in mind a possibility of wide-range charge 
doping in GdBaC^Os+z, which affects both the valence 
state of Co ions and the magnetic interactions, structural 
ordering phenomena, and a possible nanoscopic phase 
separation, in addition to the spin-state degree of free- 
dom, one might expect the magnetic phase diagram in 
this system to be extremely complicated or even messy. 
In reality, the magnetic behavior of GdBaC^Os+z turns 
out to be indeed rich, but still following fairly simple 
empirical rules. As will be discussed below, the possi- 
ble spin arrangements in GdBaCo 2 05 +2; appear to be 
dramatically simplified by an exceptionally strong spin 
anisotropy, which essentially pins the spin orientation^ 7 - 
This fortunate feature turns GBCO into a good model 
system, yet it also emphasizes the necessity to perform 
all the magnetic measurements on single crystals only. 

Figure 15 illustrates the method used to determine the 
anisotropic magnetization contribution coming from the 
cobalt ions. For every oxygen composition, we measured 
the magnetization M(T) of GdBaCo205 +2; crystals with 
the magnetic field H — 0.01;0.1;1;7 T applied along or 
transverse to the C0O2 planes. From the overall magne- 
tization we then subtracted an isotropic contribution of 
Gd 3+ ions (shown by dashed lines in Fig. 15), assum- 



ing their ideal paramagnetic behavior with spin S = 7/2; 
the latter appeared to be a good approximation since we 
found no sign of Gd 3+ spin ordering down to T = 1.7 K 
in any of the crystals. Hereafter, all the magnetization 
data will be presented after subtracting the contribution 
of Gd 3+ ions. 

When the magnetic field is applied along the 
C0O2 planes (H || ab), their magnetization exhibits 
several doping-dependent anomalies, which include a 
ferromagnetic-like behavior at moderate temperatures 
(Fig. 15). In sharp contrast, the c-axis magnetization 
(H || c) appears to be always smaller and virtually fea- 
tureless throughout the whole doping range, except for a 
weak step-like feature associated with the metal-insulator 
transition at T w 360 K for x ~ 0.5, which is observed for 
any field orientation. This persistent anisotropy clearly 
indicates that regardless of the oxygen content, the cobalt 
spins in GdBaCo205+a; are strongly confined to the C0O2 
planes. As can be seen in Fig. 15, even when cobalt ions 
develop a fairly strong ferromagnetic moment along the 
ab plane (~ 1 /xb per formula unit at T 200 — 260 K), 
the 7 T field appears to be far too weak to overcome the 
spin anisotropy and to turn this ferromagnetic moment 
towards the c-axis. In the following, we will focus on the 
in-plane magnetic behavior, since the c-axis magnetiza- 
tion is governed simply by a weak field-induced rotation 
of Co moments out of the ab plane, that is, by a compe- 
tition between the Zeeman energy and a single-ion spin 
anisotropy. 

Another general feature of the magnetization in 
GdBaCo205+a; is a pronounced thermo-magnetic irre- 
versibility. Regardless of doping, the data taken upon 
cooling the sample in magnetic field (FC) differed from 
those obtained on heating after the sample was cooled in 
zero field (ZFC), as illustrated in Fig. 16. As the mag- 
netic field increases, the ZFC curve approaches the FC 
one, yet the difference remains well discernible even at 
H = 7 T. 

Quite often such magnetic irreversibility observed in 
transition-metal oxides is attributed, sometimes without 
due care, to the formation of a spin-glass state. It should 
be noted, however, that the idea of a spin glass implies 
that no long-range (or intermediate-range) order is de- 
veloped in the spin system^ We have found that in the 
case of GdBaCo205+ 2; , the magnetic irreversibility shows 
up only below the onset of a ferromagnetic-like behavior 
(Fig. 16) and thus is most likely associated with a con- 
ventional ferromagnetic domain structure or with a meta- 
magnetic transition. It should not be surprising also if 
the FM domain structure in GBCO appears to be very ro- 
bust: one may just imagine how difficult it should be for a 
FM domain composed of Ising-like spinsi 7 - to change the 
direction of magnetization. Whatever the doping level 
in GdBaCo205+ K was, we never observed a conventional 
spin-glass behavior, and thus we suggest that all mag- 
netic irreversibilities in this system are related to the spin 
rearrangement within a spin-ordered state. 
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FIG. 16: Irreversible behavior of the magnetization in 
GdBaC^Os+rr with x = 0.33 (the contribution of Gd 3+ ions 
has been subtracted). Measurements were performed either 
on cooling the crystal in magnetic field (FC) or on heating 
after the crystal was cooled down to 2 K in zero magnetic 
field (ZFC). The magnetic field H = 0.01 T (a) or H = 0.1 T 
(b) was applied along the ab plane. 



2. Evolution with doping 

At x = 0.0, GdBaCc^Os+z exhibits an antiferro- 
magnetic behavior as was previously reported for other 
RBaCo205.o compounds studied by magnetization mea- 
surements and neutron powder diffraction)2i2i2£ the lat- 
ter identified a G-type AF ordering for spins of the cobalt 
ions [labeled as AF(1) hereafter]. It becomes increasingly 
more difficult to determine the magnetization component 
coming from the Co ions as x approaches zero, because 
the contribution from the Gd ions becomes dominant in 
the temperature range of interest; therefore, it is use- 
ful to look at the M(T) behavior of YBaCo 2 5+x - the 
only known isostructural compound with non-magnetic 
R ions"^ (Fig. 17). The magnetization of \ BaCk^O^o- 
being in rough agreement with the evaluated contribution 
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FIG. 17: Magnetization of an YBaCo205+a; ceramic measured 
at H = 0.1 T. 



of cobalt ions to the magnetization of GdBaCo 2 Os.o, is 
quite small and shows two anomalies at w 340—350 K and 
ss 220 K, which were attributed by Vogt et al. to the on- 
set of an AF spin order and a unidirectional charge order, 
respectively^ It is worth noting that the same types of or- 
dering at virtually the same temperatures were also found 
in HoBaCo 2 5 .o (Ref. 0) and NdBaCo 2 O 5 . (Ref. IH, 
which gives evidence that this spin and charge ordering 
is governed solely by the physics of Co0 2 planes. 

When the antiferromagnetic GdBaCo 2 05.o is doped 
with oxygen, it immediately causes an increase of mag- 
netization at temperatures below w 300 — 320 K [Fig. 
18(a)]; this looks like a new ferromagnetic-like compo- 
nent emerging and growing in a rough proportion to the 
oxygen concentration up to x ~ 0.3. To our surprise, the 
onset temperature of this FM behavior remains constant, 
ignoring the growth of the overall FM moment. Appar- 
ently, such behavior would certainly indicate a formation 
of a distinct chemical phase, if the structural data did not 
demonstrate that all the crystals were perfectly homoge- 
neous on a macroscopic scale (down to the scale resolved 
by conventional X-ray technique) . 

One might speculate that the doped oxygen induces 
canting of antiferromagnetically ordered cobalt spins 
which then mimic the ferromagnetic behavior, as hap- 
pens in many other canted antiferromagnets, such as 
La 2 Cu04^, where a slight spin rotation brings about 
a weak ferromagnetism. However, a rather large mag- 
nitude of FM moment, exceeding, for example, in the 
x = 0.33 sample ~ 0.3 ^b/Co at H = 7 T, is apparently 
inconsistent with such weak spin canting. 

The remaining possibility is that the ferromagnetic re- 
sponse originates from nanoscopic FM droplets imbed- 
ded in the AF(1) matrix; they may be associated with 
adjacent oxygen clusters in GdO x planes or just with 
nanoscopic regions of Co0 2 planes enriched with charge 
carriers. These FM clusters should be at least several 
unit cells in size to survive thermal fluctuations and to 
produce more than just a paramagnetic response, yet 
small enough to avoid being detected as a macroscopic 
phase. Indeed, this hypothesis finds support from the 
field dependences of the magnetization measured in mod- 
erately doped GdBaCo 2 05 +K (0 < x < 0.3), which ex- 
hibit a curious combination of ferromagnetic features (e. 
g., thermomagnetic irreversibility) with superparamag- 
netic ones. For example, in the x = 0.165 crystal, the 
M(H) curves appear to be almost perfectly linear, that 
is, M/H is independent of H , without any sign of satu- 
ration up to 7 T [Fig. 19(a)]. In fact, this is a canonical 
behavior of a superparamagnet - a paramagnetic sys- 
tem, where the role of individual spins is played by the 
moments of FM clusters. A slow saturation appears in 
the M(H) curves at higher oxygen concentrations (as an 
example, Fig. 19(b) shows M/H data for x = 0.330), in- 
dicating that FM domains grow and probably also start 
interacting with each other. 

Interestingly, the FM ordering introduced by oxygen 
seems to be unstable at low temperatures and the mag- 
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FIG. 18: Evolution of the magnetization behavior in GdBaCo205+a; with doping (the contribution of Gd 3+ ions has been 
subtracted). Selected curves are shown for three distinct oxygen-concentration regions: < x < 0.45 (a), 0.45 < x < 0.55 (b), 
and x > 0.55 (c). The M(T) data are taken on cooling in a magnetic field of 0.1 T applied along the ab plane. The inset in (b) 
shows an expanded view of the metal-insulator transition being also a spin-state transition. 



netization increase starting below 300-320 K smoothly 
gives way to an AF-like behavior as the temperature is 
further reduced [Figs. 18(a) and 19(b)]. As can be seen 
in Fig. 19(b), a strong magnetic field can affect this FM- 
AF competition, suppressing the low-temperature mag- 
netization drop: at H = 0.1 T, M diminishes by more 
than 5 times with decreasing temperature from 200 K to 
100 K, while at H = 7 T it remains almost constant. 

The sequence of magnetic transitions becomes most 
clear in the oxygen-concentration range near x ~ 0.50 
[Fig. 18(b)]. For these compositions, a net ferromagnetic 
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FIG. 19: M/H in GdBaCo20s+a: measured at several mag- 
netic fields in a range of 0.01 — 7 T for two oxygen concentra- 
tions, x — 0.165 (a) and x = 0.330 (b). The data are taken 
on cooling in a magnetic field applied along the ab plane. 



moment appears in the cobalt sublattice below « 300 K, 
and then suddenly vanishes upon further decreasing tem- 
perature (at ~ 260 K for x — 0.500), indicating successive 
PM-FM-AF transitions. The FM state shows up only 
in a narrow temperature window, which has the smallest 
width of less than 40 K for x — 0.500. Whenever the oxy- 
gen concentration deviates from the parent composition 
x = 0.500, be it on the electron-doped or the hole-doped 
side, the FM phase becomes more stable and the FM-AF 
transition shifts to lower temperatures [Fig. 18(b)]. Also, 
the low-temperature magnetization in both cases shows 
notably higher values, indicating a deviation from a pure 
AF spin order that is realized in GdBaCo205.5o- 

Another important feature of the magnetization be- 
havior in the range of 0.45 < x < 0.55 is a step- 
like change of both the magnetization and its slope at 
T w 360 K [inset in Fig. 18(b)], exactly at the tempera- 
ture of the metal-insulator transition observed in our re- 
sistivity and thermoelectric-power measurements (Figs. 
6 and 11). This change takes place in the paramagnetic 
region and is caused by the spin-state transition of Co 3+ 
ions as has been concluded based on the Curie- Weiss fit- 
ting of the PM susceptibility and structural studies of 
GdBaCo 2 5+2; and isostructural compoundsi 12 i 16 i 17 i 23 i 30 
Similar to the behavior observed in resistivity (Fig. 6), 
any deviation in stoichiometry from x — 0.50 blurs 
the spin-state/metal- insulator transition and shifts it to- 
wards lower temperatures [inset in Fig. 18(b)]. 

As the oxygen content in GdBaC^Os+a, exceeds « 
0.55, its magnetic behavior notably changes. In this mod- 
erately hole-doped region, only crystals with x ~ 0.70 
give an impression of being homogeneous [Fig. 18(c)]. 
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FIG. 20: An example of the phase separation in a 
GdBaCo2 05+! crystal with x = 0.59, which is manifested 
in the magnetization behavior as the presence of two differ- 
ent magnetic phases typical for neighboring compositions of 
x ps 0.5 and x » 0.7. Magnetization data for crystals with 
x — 0.500 and x = 0.700 are also shown for comparison. 
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FIG. 21: Magnetization of a detwinned GdBaCo205.so crystal 
in a linear (a) and logarithmic scale (b); measurements are 
done in H = 0.1 T applied along one of the crystal axes (the 
contribution of Gd 3+ ions has been subtracted). 



3. Detwinned GdBaCo2 O5.50 crystals 



For this composition, the FM order develops at much 
lower temperatures than for x < 0.55, namely, below 
150 K, suggesting a different spin arrangement of Co 3+ 
and Co 4+ ions; the FM ordering however is still followed 
by a sharp FM-AF transition at ~ 100 K [Fig. 18(c)]. 
For other oxygen concentrations, the M(T) curves clearly 
demonstrate a superposition of two magnetic phases, cor- 
responding to 1 w 0.5 and x w 0.7 (Fig. 20). What 
distinguishes this composition region from the electron- 
doped samples (x < 0.45) is that here the phase sep- 
aration (taking place in macroscopically homogeneous 
crystals^) is also well seen in the structural data, see Sec. 
IIII Al which implies a fairly large size of magnetic and 
crystallographic domains. Since the Coulomb interaction 
would prevent a formation of large charged domains, the 
phase separation should also involve a macroscopic redis- 
tribution of oxygen ions over the crystal. One might fur- 
ther speculate that the homogeneous state with x » 0.7 
corresponds to an ordered "Ortho-Ill" phase where two 
filled oxygen chains in GdO^ layers alter with one empty 
chain (O-O-X-O-O-X), in contrast to the x = 0.5 phase 
exhibiting a one-to-one alteration (O-X-O-X). It would 
be interesting to search for the superlattice peaks cor- 
responding to this putative "Ortho-Ill" structure using 
high- intensity x-ray measurements. Given a very high 
oxygen mobility, the oxygen ions could easily form or- 
dered domains when crystals were cooled down after an- 
nealing. More structural studies are however necessary 
to build a conclusive picture. 



As it becomes clear from the data presented above, the 
magnetic properties of the parent GdBaCo205.5o compo- 
sition are most important for understanding the overall 
magnetic behavior of GdBaCo 2 05 +2: . Since the crystal 
structure at x = 0.50 is orthorhombic (see Sec. IIII A|) . we 
had to detwin crystals in order to obtain single-domain 
samples for an accurate and detailed study. The magne- 
tization of a detwinned GdBaCo2 05.so single crystal (a 
fraction of misoriented domains ~ 4%) measured along 
the a, b, and c axes reveals a remarkable anisotropy of 
the spin system [Fig. 21(a)]: In the FM state that shows 
up in a narrow temperature window below 300 K, the net 
FM moment appears only along one of the orthorhom- 
bic axes, namely, along the a axis. This behavior sug- 
gests that the cobalt spins in GdBaCo205.so are not only 
strongly confined to the C0O2 planes (Sec. IIII D 1|) . but 
are also pinned to one of the in-plane directions; in other 
words, the spins system appears to be Ising-like. Note 
that a small magnetization along the b axis in Fig. 21 
comes mostly from residual misoriented domains. 

The logarithmic plot in Fig. 21(b) provides additional 
details. As we have already discussed, the c-axis magne- 
tization in GdBaCo205.5o is almost featureless similar to 
all other compositions (a small hump in the FM region 
is several hundred times smaller than for H || a; more- 
over, it may come in part from imperfect alignment of 
the crystal in the magnetometer). The behavior of the 
in-plane magnetization is more interesting: In the FM 
region, M is much larger for H || a than for H || b, but 
the situation abruptly turns over as the crystal enters 
the AF region. Given that the transverse susceptibility 
of an antiferromagnet exceeds the longitudinal one, we 
can conclude that the cobalt spins keep their spin easy 
axis (|| a) upon the FM-AF transition. In other words, 
the cobalt spins are aligned along the a axis in the FM 
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FIG. 22: (a) a-axis magnetization of GdBaCo2 0s.5o measured in different magnetic fields, (b) AF-FM phase boundary 
determined from the M(T) and M(H) data (open circles) and magnetoresistance measurements (solid circles), (c) Isothermal 
magnetization M{H) measured upon increasing magnetic field H || a at several temperatures; the AF-FM transition exhibits a 
weak hysteresis (not shown). 



state and they keep being aligned along the a axis in the 
AF state; what happens upon the FM-AF transition is 
that the spins in one of sublattices just flip by 180°. 

The balance of FM and AF ordering in GdBaC^Os.so 
turns out to be quite delicate, so that it can be easily 
affected by magnetic fields applied along the spin-easy a 
axis, which stabilize the FM state and shift the FM-AF 
transition to lower temperatures (Fig. 22). We find that 
the magnetic field H || a required to overcome the AF 
coupling grows roughly linearly with cooling, from zero 
at T » 260 K up to ~ 20 T at T = 0. It is worth not- 
ing that the AF-FM transition remains very sharp even 
in the temperature range close to 260 K, where rather 
weak fields (~ 1.5 T at 240 K, i.e. ^ B H < kT) are ca- 
pable of recovering the FM order [Figs. 22(a) and 22(c)], 
which shows that thermal fluctuations are irrelevant here. 
This behavior clearly indicates that the observed AF-FM 
switching is a metamagnetic transition, that is, it occurs 
within the ordered spin state and is governed by the rel- 
ative reorientation of weakly-coupled spin sublattices. 

Whatever the temperature, the 7-T field H || 6 or 
H || c appears to be too week to compete with the spin 
anisotropy; it neither can turn the FM moment towards 
the b or c axis in the FM region, nor can it induce the AF- 
FM transition at T < 260 K. The only impact of H || be is 
therefore to cause a partial tilting of the cobalt spins from 
their easy a axis, thus giving linear M(H) curves in both 
FM and AF regions (Fig. 23). The spin anisotropy be- 
tween the c axis and the ab plane is several times stronger 
than the in-plane anisotropy; consequently, the slope of 
the M(H) curves for H || c is several times smaller than 
that shown in Fig. 23. 

It is interesting to examine the magnitude of the mo- 
ment that shows up in the FM state, which can be easily 



done using the M{H) data in Fig. 22(c). The FM mo- 
ment turns out to exceed 0.6 /^b/Co at T = 200 K, and a 
rough extrapolation to T = suggests a saturated mag- 
netic moment of ~ 1 /is/Co (Fig. 24). It is worth noting 
that poly crystalline RBaC^Os.s samples were reported 
to demonstrate noticeably smaller FM moments ^ii&2Si 
This apparent discrepancy originates from the Ising spin 
anisotropy that prevents moments from being seen along 
the b and c detwinned single crystal is clearly 

necessary to make the true FM moment of ss 1 /xb/Co 
visible. 

While this value is apparently too large to be ac- 
counted for by any spin-canting picture, it agrees well 
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FIG. 23: Isothermal magnetization of GdBaCo2 0s.5o for 
H || b. Dashed lines indicate the experimental data, and sym- 
bols show M(H) after subtracting the contribution from mis- 
oriented domains (4% of the total amount). 
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FIG. 24: Net magnetic moment along the a axis measured at 
H = 7 T (solid squares), and the extrapolated FM moment 
at H = (open squares). Dashed lines are a guide to the eye. 



with the expectation for a simple FM order, if the Co 3+ 
ions in this temperature range possess a 1:1 mixture of 
low-spin S = and intermediate-spin S = 1 states. In 
fact, the same conclusion on the spin states of Co 3+ ions 
has been also reached based on the Curie- Weiss fitting 
of the PM susceptibilitjiiSiiSiii in the temperature range 
300-360 K (Figs. 18 and 22), as well as based on the 
structural data^S 
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FIG. 25: (a) Temperature dependence of p c measured at H = 
and 14 T applied along the ab plane. The dashed lines show 
simple activation fits p c oc exp(A/T) for both the AF and 
FM states, and Tn indicates the temperature of the FM-AF 
transition that is shifted to lower temperatures by the applied 
field, (b) A linear-scale view of the high temperature region, 
illustrating a step-like increase of p c at Tn. 



E. Magnetoresistance 

In compounds with competing magnetic orders, a mag- 
netic field that favors one kind of ordering often causes 
also a large magnetoresistance; GdBaCo20s +2: is no ex- 
ception. The charge transport in this system turns out 
to be very sensitive to both the FM and AF ordering, 
and magnetic fields readily induce a giant magnetoresis- 
tance by affecting the subtle AF-FM balancei^S*^ As an 
example, Fig. 25 shows the c-axis resistivity of a twinned 
GdBaCo205.5o crystal measured at H = and 14 T. In 
zero field, the FM-AF transition at « 260 K brings about 
a step-like increase of the resistivity p c [Fig. 25(b)]; a 
similar, albeit smaller, step is observed in p ao as well. A 
14-T field applied along the ab plane shifts the magnetic 
transition towards lower temperatures and wipes out the 
resistivity increase, thus causing the resistivity to drop 
by up to several times. 

As can be seen in Fig. 25(a), p c grows roughly expo- 
nentially upon cooling below ~ 230 K regardless of the 
applied field; the charge carriers, therefore, have differ- 
ent activation energies in the FM and AF states, and the 
MR originates from the reduction in this activation en- 
ergy. As soon as the magnetic field becomes insufficient 
to maintain the FM order, the system switches into the 
AF state [at T = Tn(14T)], and the resistivity jumps to 
its zero-field value [Fig. 25(a)]. 

Unfortunately, we did not succeed in detwinning the 
crystals with already prepared electrical contacts, and 
thus the measurements were carried out on twinned crys- 
tals, which does not allow us to analyze the MR quanti- 
tatively. The role of twins is less critical for p c (T), since 
we find the orthorhombic domains to always go through 
the whole crystal from one face to another. Consequently, 
the measuring current does not cross domain boundaries, 
and the observed magnetoresistance Ap c / p c is just re- 
duced from its true value by some factor, since not all 
domains are affected by the magnetic field. Upon mea- 
suring p a b, however, the current flowing in the ab plane 
has to pass through both kinds of orthorhombic domains 
which form a striped structure. The measured Ap ao / p ao 
thus depends not only on the ratio of domains, but also 
on the yet unknown anisotropy p a / Pb- 

In spite of this complication, the qualitative behavior 
of the MR remains clear. At T < T N w 260 K, the in- 
plane magnetic field induces an abrupt decrease in both 
p c and pab (Fig. 26), which occurs at exactly the same 
field as the AF-FM transition in magnetization, leaving 
no doubts about its origin. At higher temperatures, the 
resistivity changes gradually, again resembling the mag- 
netization behavior in Fig. 22(c). In the latter case, 
the MR seems to originate from the field-induced stabi- 
lization of the FM order; the magnetic field suppresses 
critical FM fluctuations, thus facilitating the charge mo- 
tion. 

The temperature dependences of Ap ao / p ao and Ap c /p c 
in Fig. 26(b) nicely illustrate the correlation of the MR 
with the magnetic behavior. Upon cooling, we first ob- 
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FIG. 26: Magnetoresistance of twinned GdBaCc^Os.s crystals: field dependences of Ap c /p c (a) and Ap a t>/pab (c) measured at 
several temperatures for H || ab. For clarity, in panel (c) the hysteresis is shown only for T = 120 K; for other temperatures, 
the data are given for sweeping the field one-way from -14 T to +14 T. (b) Temperature dependence of Ap c /p c and Ap a b/p a b 
measured at H = 14 T applied along the ab plane. 



serve a clear peak at ~ 278 K in the MR, which can 
be associated with the FM transition (the peak position 
coincides with the peak in dM/dT in Fig. 21). Ap- 
parently, the FM fluctuations that are strongest near 
Tq frustrate the charge motion, providing contribution 
of about 20% to resistivity, and this contribution is re- 
moved when a strong magnetic field is applied. It is worth 
noting, that this negative MR peak near Tq closely re- 
sembles the MR behavior of such ferromagnetic oxides 
as cubic Lao.sBao.sCoOa (Ref. 0) and Lai-^Sr^MnOa 
(Ref.[l|), where the resistivity change is usually explained 
by the double-exchange mechanism^ Upon further cool- 
ing, the MR abruptly increases at T = Tn- The ori- 
gin of this MR is obvious: Without field, the AF or- 
dering enhances the resistivity, while the applied mag- 
netic field prevents establishing the AF order. As can be 
seen in Fig. 26(b), the c-axis MR quickly gains strength 
with cooling. At 150 K, p c drops by three times at 
14 T, and the ratio p c (0,T)/ p c (H,T) would keep grow- 
ing to much larger values with lowering temperature if 
fields exceeding H C (T = 0) were applied^ The in-plane 
MR Ap a b/pab seems to be by several times smaller; the 
magnitude and the temperature dependence of Ap a b/p a i, 
should, however, be taken with a grain of salt, since, as 
we mentioned above, a detwinned crystal should be mea- 
sured to determine the true MR values. 

The MR anisotropy with respect to the magnetic- 
field direction can provide interesting information on the 
cobalt-spin anisotropy. In general, when an Ising-like an- 
tiferromagnet is subjected to an increasingly strong mag- 
netic field, it eventually turns into a FM state, whatever 
the field direction is. If the field is applied along the spin- 
easy a axis, the spin-flip transition is abrupt, while for the 
transverse direction, the spins rotate gradually and much 



higher fields are required to align them. In the latter 
case, the rotation angle is roughly sin a ~ gp^H/2 J a _ c , 
where J a _ c is the spin anisotropy energy. It is reasonable 
to expect that the MR in GdBaCo205 +:c , being roughly 
proportional to the magnetization, will change with the 
field as cx gp^H/2 J a _ c , tending to reach its saturation 
value at H ~ 2J a _ c / gpB- Thus, the anisotropy of mag- 
netoresistance in addition to that of magnetization may 
be used to probe how strongly the spins are coupled with 
the crystal lattice. 

An experimental study of GdBaC^Os.so crystals re- 
veals that the MR anisotropy is surprisingly large and 
the 14-T field H || c can do nothing comparable to the 
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FIG. 27: Dependence of the magnetoresistance on the field 
direction: Ap c /p c (left) and Ap a i,/p a t (right) measured at 
T = 240 K with the magnetic field applied along or transverse 
to the ab plane. 
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MR caused by a much less field H || ab (Fig. 27). Indeed, 
the MR is barely seen for H || c, especially Ap ab /p ab [Fig. 
27(b)], which is as small as 0.14% at 14 T, while for the 
in-plane field it readily reaches a two orders of magni- 
tude larger value. This indicates that magnetic fields in 
the 100- T range would be necessary to rotate the cobalt 
spins in GdBaCo 2 05.5o from the a to the c axis, over- 
coming the remarkable spin anisotropy. It is worth not- 
ing, that such spin-anisotropy energy, roughly estimated 
to be of the order of 10 meV/Co (fi B H), is extremely 
large compared with any known magnetic materialise 
and is close to the giant magnetic anisotropy observed 
in nanoclusters,— This huge spin anisotropy can also 
account for a considerable anisotropy in susceptibility, 
Xab/Xc > 1) that survives even in the paramagnetic state 
up to 400 K. 



F. Phase diagram of GdBaCo2 05+a; 

An empirical phase diagram of GdBaCo 2 05+ 2; based 
on the obtained structural, transport, and magnetic 
data is sketched in Fig. 28 as a function of the 
oxygen concentration. At the lowest oxygen content, 
x = 0, GdBaCo205+ x is an antiferromagnetic insulator 
[AFI(l)]; its low-temperature phase should most likely 
correspond to a charge-ordered G-type antiferromagnet, 
by analogy with other isostructural compounds studied 
using neutron diffractioni&2i5i 

An increase in the oxygen content x is found to imme- 
diately result in the formation of isolated ferromagnetic 
clusters imbedded in the antiferromagnetic [AFI(l)] ma- 
trix. These clusters demonstrate a paramagnetic behav- 
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FIG. 28: Phase diagram of GdBaC^Os+o;, including re- 
gions of a paramagnetic metal (PMM), paramagnetic insu- 
lator (PMI), ferromagnetic metal (FMM), ferromagnetic in- 
sulator (FMI), and antiferromagnetic insulator (AFI). 



ior above Tq = 250 — 300 K (marked by open triangles 
in Fig. 28), and a curious combination of superpara- 
magnetic and ferromagnetic features below 7c; namely, 
the magnetization does not saturate in magnetic fields 
up to 7 T (as in canonical superparamagnets), yet the 
FM moments exhibit an extremely strong coupling with 
the crystal lattice, which is manifested in the magnetic 
anisotropy and thermo-magnetic irreversibility. 

With further increasing oxygen content in 
GdBaCo 2 5+a; (up to x = 0.45), the FM clusters 
smoothly grow in size, reaching percolation, and develop 
at low temperatures an intrinsic instability towards a 
new type of AF ordering [AFI(2)], which is the ground 
state of the parent GdBaCo20s.5o compound^ The 
AFI(2) order emerges from the FM order in clusters and 
differs essentially from the AFI(l) state that is realized 
at x — 0; for instance, the AFI(2) ordering can be easily 
suppressed by a magnetic field, while AFI(l) is quite 
robust. One can consider GdBaCo205-|_ :r over a wide 
composition region of < x < 0.45 as being composed 
of nanoscopic phases mixed together: At high tempera- 
tures, the system behaves as a paramagnetic insulator; 
then, below Tq, as a mixture of FM and AF insulating 
phases; and finally at low temperatures, it evolves into a 
mixture of two distinct AF phases [AFI(l) and AFI(2)]. 
The resistivity in this composition range shows a robust 
insulating behavior without any noticeable change with 
increasing x. The room-temperature crystallographic 
structure retains a macroscopically tetragonal symmetry 
with smoothly changing lattice parameters. 

The region with the richest behavior, 0.45 < x < 0.55, 
occupies the center of the phase diagram around the 
"parent" composition x = 0.50. In this doping range, 
the crystal structure is macroscopically orthorhombic up 
to rather high temperatures due to the ordering of oxygen 
into alternating full and empty chains that run along the 
a axis (at T — 260°C, we still observed a twin structure 
related to the oxygen ordering) . Owing to this structural 
order, the broad magnetic and transport features emerg- 
ing already at x = 0.3 — 0.4 come in focus here, and 
GdBaCo205+2; demonstrates a series of sharp phase tran- 
sitions upon cooling: first, from a paramagnetic metal 
(PMM) to a paramagnetic insulator (PMI) at Tmit (open 
squares in Fig. 28), which is also accompanied with a 
spin-state transition, then to a ferromagnetic insulator 
(FMI) at Tc, and, finally, to an antiferromagnetic insu- 
lator [AFI(2)] at Tfm-AF (open circles). 

The "parent" composition GdBaCo 2 05.5o is a metal at 
high temperatures but turns into a semiconductor upon 
cooling below the metal-insulator transition. In fact, the 
x = 0.50 composition represents a kind of borderline that 
divides the central region of the phase diagram into two 
roughly symmetric parts which correspond to electron- 
doped and hole-doped semiconducting states below Tmit- 
All transport properties show anomalies in their doping 
dependences at x = 0.50, indicating that the density of 
doped carries smoothly goes to zero with approaching 
this point from either side. As can be seen in Fig. 28, all 
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the transition temperatures also exhibit maximum val- 
ues exactly at x = 0.5. The only property demonstrating 
a striking asymmetry with respect to x — 0.5 is the re- 
sistivity which is much smaller for the hole-doped side 
(x>0.5). 

When the oxygen concentration exceeds x — 0.55, 
GdBaCo 2 05 +2; develops a new phase possessing consid- 
erably lower temperatures of the ferromagnetic transition 
(T c < 150 K) and the FM-AF transition (T F m-AF < 100 
K). Until x reaches « 0.7, this new phase is mixed (on 
a mesoscopic scale) with the x ~ 0.5 phase, and only at 
x > 0.68—0.70 the system again recovers its homogeneity. 
Surprisingly, in the x w 0.7 region, the series of succes- 
sive PM-FM-AF transitions, albeit happening at reduced 
temperatures, still remains remarkably similar to that at 
x ks 0.5, in spite of very different transport properties 
and different spin states of cobalt ions in these phases. 

In general, GdBaCo205+ x seems to have just a few 
stable phases such as x = 0, x » 0.5, and x « 0.7, while 
intermediate compositions always tend to phase separate 
on a nanoscopic or mesoscopic scale. 

Although the compositions with x > 0.77 are hard to 
achieve in GdBaCo205+ x , we can extrapolate the phase 
boundaries to higher doping, as shown in Fig. 28, us- 
ing an analogy with the cubic Lao.sBao.sCoOs-^. Our 
measurements have shown that the latter compound is 
a ferromagnet with Tq « 200 K (only weakly dependent 
on S) with no sign of reentrant AF behavior at low tem- 
peratures. 



IV. DISCUSSION 

It is natural to start discussing the GdBaCo205+ x 
compound beginning with the "parent" x — 0.50 com- 
position, whose properties constitute a basis for under- 
standing the transport and magnetic behavior of this lay- 
ered cobalt oxide over the entire doping range. 

A. Magnetic structure of GdBaCo^Os.so 

By now, there were several attempts to eluci- 
date the magnetically ordered states in isostructural 
RBaCo 2 5+2; (R = Y, Nd, Tb, Ho) compounds us- 
ing neutron scattering (for GdBaCo205 +a; , such a study 
is precluded by a large absorption of neutrons by 
Gd)i 8 i 9 i 51 i 58 i 59 i 60 i 61 In the case of samples with the low- 
est oxygen concentration x — 0, whose preparation is 
quite straightforward, the neutron diffraction has indeed 
provided reproducible data on the charge and spin or- 
dering, being in good agreement with the magnetiza- 
tion and transport measurements In contrast, much 
more complicated and sample-dependent diffraction pat- 
terns were obtained for compositions near x « 0.5, bring- 
ing about controversial models of the magnetic order in 
this region of the phase diagrami^SiffiiSi This apparent 
controversy, however, is in fact not very surprising given 



the remarkable sensitivity of the magnetic and structural 
properties of RBaCo205+ :c to even a slight modification 
in the oxygen concentration, that was revealed in the 
present study; note also a strong tendency to phase sepa- 
ration, which should inevitably take place in all the sam- 
ples except for those located in a few narrow composi- 
tion regions (see Fig. 28). As we have shown here, the 
actual composition of so-called "as-grown" and "oxygen- 
annealed" samples usually used in most of the previous 
studies may deviate a lot from the required x = 0.50, and 
a rather delicate technique should be developed to tune 
the oxygen content precisely to this peculiar point. 

To the best of our knowledge, none of the models, that 
have been suggested for the x — 0.50 composition based 
on the neutron-scattering data, appears to be capable of 
explaining the entire experimentally-observed magnetic 
behavior. For example, a G-type AF order was proposed 
by Fauth et alj& for temperatures where a rather large 
FM moment shows up in our magnetization measure- 
ments. On the other hand, the model suggested by Soda 
et a/.|S2. being quite close to that one reported recently by 
uSfli indeed captures such gross features as the FM and 
AF states and a possibility of easy and abrupt switching 
between these two; nevertheless, it still implies a kind of 
spiral spin order for the AF state, which is clearly incon- 
sistent with the Ising-like magnetic behavior observed in 
the magnetization of detwinned crystals^ 

In view of the above mentioned problems with di- 
rect methods, we should try to reconstruct the magnetic 
structure of GdBaCo205.5o relying mostly on the mag- 
netization data obtained on single crystals with precisely 
tuned stoichiometry. Fortunately, this task is simplified 
by the discovered strong uniaxial anisotropy of cobalt 
spins which dramatically narrows down the range of pos- 
sible spin arrangements. Let us first summarize the ro- 
bust facts established for GdBaCo20s.5o as follows: 

(i) The magnetic moments of Co 3+ ions in 
GdBaCo205.5o exhibit an Ising-like behavior, being ex- 
ceptionally strongly confined to the C0O2 (ab) planes 
(with the anisotropy energy of the order of 10 meV per 
cobalt ion) and to a lesser extend, though still strongly, to 
the orthorhombic a axis. Note also that the confinement 
of cobalt spins to the ab plane appears to be a generic 
feature of GdBaCo205 +a; regardless of the oxygen con- 
tent and thus regardless of the valence and spin states 
of cobalt ions fSec. IIII D 1|) . One might wonder whether 
there can be another interpretation of the observed mag- 
netic anisotropy: In magnetically ordered states, for ex- 
ample, the magnetization may indeed be larger for a di- 
rection transverse to the spin-easy axis, which is usually 
the case with antiferromagnets and canted antiferromag- 
nets. However, a possibility of the spin-easy axis to be 
parallel to the c axis is immediately ruled out by a con- 
siderable anisotropy Xab > Xc that survives up to high 
temperatures in the paramagnetic state. In turn, by an- 
alyzing various canted AF configurations within the ab 
plane, regardless of whether they are allowed by symme- 
try or not, we find that non of them can account for the 
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observed large FM moment, which would require a large 
canting angle of cobalt spins, and at the same time for 
a strong in-plane anisotropy, where the FM moment ap- 
pears only along the a axis (Fig. 21). For example, the 
spin structure suggested in Ref . |59j with a canting angle 
of ~ 7r/4 would inevitably result in a magnetic behavior 
to be essentially isotropic within the ah plane. 

(ii) The spontaneous FM state below Tq and the 
FM state induced by a magnetic field H || a at T < 
Tfm-af represent the same phase with a net FM moment 
smoothly growing upon cooling, tending tow 1 Mb/ Co at 
T = (Fig. 24). After the magnetic field stabilizes the 
FM order, no anomaly can be seen in the magnetization 
or resistivity at the zero-field FM-AF transition temper- 
ature. 

(hi) The FM-AF switching in GdBaCo 2 05. 5 o is a meta- 
magnetic transition: It is induced by a relative reorienta- 
tion of weakly coupled spin sublattices, while within each 
sublattice the spins are kept ordered by a much stronger 
interaction. This remarkable hierarchy of spin interac- 
tions follows from the fact that the AF-FM transition 
remains sharp even when it is induced by a very weak 
field fi B H < kT (Fig. 22). 

Based on the above listed observations (i)-(iii), one can 
quite easily sort out the spin structures that may be re- 
alized in GdBaCo2C>5.5o. The point (i) indicates that the 
cobalt spins in this compound should form a collincar 
structure aligned along the a axis; this conclusion agrees 
with all the neutron-scattering studiesj2i2i£ii£2iffi except 
for Ref. However, if we consider a uniform FM order 
at T < Tq, the observed net FM moment of w 1 /ib/Co 
[see (ii)] appears to be inconsistent with the allowed spin 
states of Co 3+ ions (0, 2/zb, and 4/xb for LS, IS, and HS 
states, respectively). Although the magnetic moments of 
Co 3+ ions could be modified from their spin-only values 
by the orbital moments, which at a first glance may re- 
solve this contradiction, a homogeneous FM state is also 
clearly inconsistent with the FM-AF switching (iii). We 
thus inevitably arrive at a conclusion that the magnetic 
structure of GdBaCo20s.5o should correspond to a fer- 
rimagnet composed of two or more antiparallel FM sub- 
lattices; note that one of the sublattices may carry zero 
moment if it involves the LS cobalt ions. Apparently, 
the net FM moment of 1 fiB /Co - a difference between 
the sublattice moments - may be accounted for if the 
Co 3+ ions constituting the two ferrimagnetic sublattices 
are either in the LS and IS states, or in the IS and HS 
states. One can construct more complex spin structures 
by increasing the number of magnetic sublattices; for in- 
stance, a zero-moment sublattice may be composed of LS 
Co 3+ ions, or antiferromagnetically ordered IS(HS) ions. 

Now we can try to map the ferrimagnetic sublattices 
onto the actual crystal structure, keeping in mind the 
following important conditions. First, the ferrimagnetic 
sublattices should be strongly coupled with each other, 
since, according to Respaud et at, even a 35-T field fails 
to induce a ferrimagnetic-ferromagnetic transition^ An- 
other explanation for this robustness is that only one of 



the sublattices actually possesses a non-zero FM moment; 
in that case, no decoupling transition should be expected. 
The second point is that besides strong FM bonds, each 
3D ferrimagnetic sublattice should also contain planes of 
very weak magnetic coupling to allow for an easy fold- 
ing of the magnetic unit cell upon switching into the AF 
state. One can see in Fig. 1, that the crystallographic 
planes responsible for weak magnetic coupling can nat- 
urally be the GdO^ plane and/or an ac plane passing 
between the pyramidal and octahedral cobalt ions. Ap- 
parently, such kind of effective decoupling of neighbor- 
ing Co 3+ ions can hardly be possible unless either of 
the two ac C0O2 planes - composed of cobalt ions in 
the pyramidal or octahedral environment - is essentially 
non-magnetic. Consequently, the other half of cobalt ions 
should provide an average moment of ~ 2 /ib/Co, which 
corresponds to the IS state or an w 1:1 ratio of LS and 
HS states. Note that the existence of a non-magnetic 
sublattice also gives the best account for the high-field 
magnetization data of Respaud et al.J& 

An important information on the magnetic states of 
Co 3+ ions in GdBaCo20s.5o has been also obtained from 
the Curie- Weiss fitting of the PM susceptibility 12 i 16 i 17 
in the temperature range 300-360 K (Figs. 18 and 22), 
which better agrees with the presence of 50% of ions in 
the IS state, while the same fitting above Tmit ~ 360 K 
points to a 1:1 mixture of IS and HS states. Although our 
magnetization experiments can hardly distinguish which 
cobalt ions alter their spin state at the spin-state transi- 
tion at Tmit and whether the IS state is realized in pyra- 
midal or octahedral positions, a structural analysis of the 
oxygen coordination by Frontera et al. suggests that it 
is the cobalt ions in octahedral positions that participate 
in the spin-state (LS-HS) transition, while the cobalt 
ions in pyramidal sites keep their IS state regardless of 
temperature^ 

We have eventually reached the most likely picture of 
the magnetic ordering in GdBaCo205.so, which is illus- 
trated in Fig. 29. The alternating filled and empty oxy- 
gen chains shown in Fig. 29(a) create two types of struc- 
tural environment, octahedral and pyramidal, for Co 3+ 
ions; the former favors the non-magnetic LS ground state, 
while the latter makes the IS state preferable. At high 
temperatures, the entropy keeps the cobalt ions in oc- 
tahedral positions in the HS state, and the compound 
behaves as a paramagnet composed of IS and HS Co 3+ 
ions (not to mention Gd 3+ ions which are always para- 
magnetic). Upon cooling below Tmit ~ 360 K, the octa- 
hedral cobalt ions cooperatively switch into the LS state, 
and thus the C0O2 planes develop a spin-state order con- 
sisting of alternating rows of Co 3+ ions in the LS and IS 
states. Consequently, the magnetic cobalt ions form 2- 
leg ladders extended along the a axis and separated from 
each other by non-magnetic C0O2 ac layers along the b 
axis and by paramagnetic GdOo.5 layers along the c axis, 
as sketched in Fig. 29(a). 

Empirically, the spin interaction in ladders formed be- 
low Tmit turns out to be ferromagnetic, and eventually 
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FIG. 29: A sketch of the magnetic structure of GdBaC02O5.sc) in the FM state (a) and possible AF states (b)-(d); Ba and Gd 
ions are omitted for clarity. At T < Tmit, the Co 3+ ions in octahedral environment are considered to have a low-spin (S = 0) 
state, while those in pyramidal environment have an intermediate (S = 1) state (the spin direction is indicated by arrows). 



the ladders establish a FM order below T c « 280 - 300 
K with spins aligned strictly along the a axis [Fig. 
29(a)]. Owing to the reduced dimensionality (quasi- 
1D/2D) of the ladders, the FM order develops quite 
gradually upon cooling, being subject to strong ther- 
mal fluctuations. The strength of the FM interaction 
in ladders, J/k-Q, can be estimated from the Curie tem- 
perature: The molecular-field theory for 3D systems^ 
would give J/fcs ~ 150 K for Tc ~ 300 K; however, for 
the quasi- 1D/2D magnetic ordering, J should be roughly 
twice as large for the same Tc to be reached** 3 - It is 
worth noting that GdBaCo205.5o is an insulator and 
thus such mechanism of the FM ordering as the double 
exchange^ which is common in metallic oxides^ is irrel- 
evant here. In insulators, the magnetic ordering is known 
to be predominantly caused by the superexchange inter- 
action, whose sign for each pair of ions can be estimated 
using Goodenough-Kanamori ru l eSl £4i§iiiSii§I To account 
for the FM spin order in ladders, we need to assume some 
kind of orbital ordering among the IS Co-ions>i£ We can 
even speculate that the superstructure reflections along 
the a axis often observed in the neutron scattering are 
brought about by such orbital ordering. 

Whether a macroscopic magnetic moment will emerge 
or not after the FM-ordered ladders are formed depends 
on their relative orientation, which can be either ferro- 
or antiferromagnetic due to the Ising-like nature of the 
spins. In a narrow temperature range below Tc, the weak 
effective interaction between ladders turns out to be fer- 
romagnetic and a net FM moment shows up; upon fur- 
ther cooling below Tfm-AF, however, it gradually changes 
sign, resulting in the AF ground state. The ladder stack- 
ing may become antiferromagnetic along the b axis [Fig. 
29(b)], inducing doubling of the magnetic unit cell along 
that axis, or along the c axis [Fig. 29(c)], or both [Fig. 
29(d)]. 

Quite naturally, the inter-ladder coupling across non- 



magnetic layers is weak, and switching from the ground- 
state AF order to the field-induced FM one can be in- 
duced by a fairly weak magnetic field H c that depends 
on temperature. To understand the role of temperature, 
one should consider the thermally-excited states: upon 
heating, a certain amount of LS Co 3+ ions in octahe- 
dral positions become Co 2+ , Co 4+ , or change their spin 
state. Whatever the case, each excited ion acquires a 
non-zero spin, providing an additional bridge between 
spin-ordered FM ladders. In this respect, it is worth not- 
ing that the AF ordering can be suppressed not only by 
magnetic fields or thermal excitations, but also by chang- 
ing the oxygen stoichiometry: any deviation from x = 0.5 
also introduces Co 2+ or Co 4+ ions and shifts the AF-FM 
phase transition to lower temperature [see Fig. 18(b) and 
Fig. 28]. 

We can conclude that the GdBaCo205.so compound 
with precisely tuned stoichiometry represents a very in- 
teresting magnetic system of weakly interacting FM lad- 
ders with Ising-like moments. 



B. Electronic structure of GdBaC^Os.so 

As follows from the magnetic and structural data of 
GdBaCo205.5o discussed above, the metal-insulator tran- 
sition at T MIT w 360 K is accompanied with a change in 
the spin-state of Co 3+ ions: Upon heating across Tmit, 
the Co 3+ ions in octahedral positions switch their spin 
state from the LS to HS state, while those in pyramidal 
positions keep their IS state at all temperatures. Fig- 
ure 30(a) shows the electronic energy levels of Co 3+ ions 
for the relevant states, namely, the LS and HS states in 
octahedral environment, and the IS state in pyramidal 
environment. This level structure is evaluated in the 
framework of the ionic model that takes into account 
the crystal-field (CF) energy splitting (10D g ) and on-site 
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FIG. 30: Schematic representation of the electronic struc- 
ture of GdBaCo2O5.50- (a) Electronic levels of a Co 3+ ion in 
the LS state (octahedral environment); IS state (pyramidal 
environment); and HS state (octahedral environment). The 
electronic band structure of C0O2 planes below the metal- 
insulator/spin-state transition originates from a superposition 
of the LS(oct.) and IS(pyr.) states (b), while above the tran- 
sition it originates from IS(pyr.) and HS(oct.) states (c). 



electrons interactions (the intra-atomic exchange energy 
Jh and the Coulomb repulsion energies U and U', defined 
following Kanamor&S) . 

In the octahedral environment (Oh symmetry), the 
crystal field splits the 3d electron energy levels of a 
Co 3+ ion into three ti g orbitals (d xy , d xz , and d yz with 
AE = —AD q ) and two e g orbitals (d x 2_ y 2 and d 3z 2_ r 2 
with AE = +6D q ). When such cobalt ion acquires 
the LS state, all its six valence electrons occupy the de- 
generate tig orbitals (three d xyxzyz spin-up and three 
dxy,xz,yz spin-down). Every electron therefore has the 
same energy, and the total electron energy per Co 3+ 
ion, which takes into account the on-site interactions and 
the crystal field splitting, can be expressed as £^ c f = 
E - 24D q - 6 J H + 3U + 12U', where E is a reference 
energy. 

In order to bring the Co 3+ ion into the IS state, one 
electron should go from a tin orbital to a higher-located 
e g one, for example, from dj y to d^ z2 _ r2 . This leads to 
an additional splitting of the one-electron energy levels 
that removes the initial degeneracy of ti g orbitals: the 
intra-atomic exchange is modified, causing every spin- 
up electron to lower its energy by Jh, while every spin- 
down electron increases the energy by Jh- Besides, the 
Coulomb matrix element U between the d xy and d^. y or- 



bitals is replaced with U' between the d xy and d^ z2 _ r2 
orbitals. Consequently, the total electron energy for the 
IS state becomes Ej£ = E - UD q - 7 J H + 2U + UU'. 
And finally, in the case of the Co 3+ HS state in the oc- 
tahedral environment, one more electron is transferred 
from the ti g to e g shell (from d@ z to d" 2 _ y2 in Fig. 30(a)), 
causing a further splitting of the energy levels. The to- 
tal electron energy for the HS state can be written as 



f hs 

^oct 



E - AD q - 10 J H + U + UU'. 



The energy levels of cobalt ions in the pyramidal envi- 
ronment (C± v symmetry) differ from the previous cases 
because the crystal field causes a different energy split- 
ting: +9.14D ? (d x 2_ y2 ), +0.86D q (d 3z 2_ r2 ), -0.86D, 
{d xy ) and -i.57D q (d xz and d yz ) for a square pyramid^ 
Consequently, the total electron energy for the IS state, 
for instance, becomes E™ r =E - 18.28 L> 9 - 7 J H + 2U + 
UU'. 

In order to evaluate which spin states of Co 3+ ions 
should be preferable, we need to compare their relative 
energies. Given that the difference between the intra- and 
inter-orbital Coulomb matrix elements is U — U' = 2Jh, 
only two parameters Jh and D q appear to determine the 
one-electron energy levels of localized Co 3+ states in the 
framework of the ionic model. Quite reasonable param- 
eters for Co 3+ ions are Jh ~ 0.5 eV and D q ~ 0.25 
eVpS Following this approach, one can easily find out 
that the IS state in octahedrons is completely unstable 
with respect to either the LS or HS state, while either of 
the latter states may win, depending on the parameters. 
However, for the pyramidal positions, the IS state ap- 
pears to be stable in a very broad range of J\i/D q (from 
0.57 to 2.74), as has been concluded by Pouchard et al.m& 

It is worth noting also that the result of this spin- 
state competition, that is, which one of the spin states 
closely located in energy will win, may well depend on 
temperature. For example, the non-magnetic LS state 
in octahedral positions experimentally appears to be the 
ground state, and therefore it dominates at low temper- 
atures. However, the spin-orbital degeneracy of the HS 
state is much larger, 2 ? and thus with increasing temper- 
ature, the HS state should sooner or later take over the 
LS state because of the entropy terms associated with the 
spin-orbital degeneracy; this mechanism is likely staying 
behind the spin-state transition observed at Tmit ~ 360 
K. 

Another important issue is the energy-level position of 
occupied e g orbitals. As shown in Fig. 30(a), for IS Co 3+ 
ions located in the pyramidal surrounding, the energy of 
d 3z 2_ r 2 occupied orbital is very close to the energy of d xy , 
d xz , and d yz orbitals, mainly due to the strong splitting 
of d 3z 2_ r 2 and d x 2_ y 2 levels by the crystal field of the Ci v 
symmetry. In contrast, for cobalt ions in the octahedral 
environment, the energy of d 3z 2_ r 2 and d x 2- y 2 occupied 
orbitals appears to be much higher even in the HS state. 

Quite naturally, the qualitative features of these sim- 
ple ionic model may be expected to be reproduced in the 
band structure. In the framework of a more sophisti- 
cated ligand-field (LF) approach, the overlap of cobalt 
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3d, 4s, 4p and oxygen 2d orbitals leads to the formation 
of bonding, antibonding, and non-bonding molecular or- 
bitals (MO), which are broadened in a solid into energy 
bands. The bonding molecular orbital, being predomi- 
nantly composed of the s and p "atomic orbitals" (AO), 
creates a wide valence band (indicated in Figs. 30(b) and 
30(c) as V.B.) similar to ordinary band insulators (semi- 
conductors). In turn, the antibonding molecular orbital, 
being also predominantly composed of s and p AO of 
cobalt and oxygen ions, creates a wide conduction band 
(indicated as C.B.). 

The interaction of cobalt d orbitals with oxygen p or- 
bitals has a different character. In a C0O6 octahedron, 
the cobalt i 2 g orbitals should be non-bonding, because 
there is no linear combination of oxygen orbitals of the 
same symmetry for them to interact. As a result, a rather 
narrow ir* band occupied by almost localized electrons 
appears at the top of the valence band. On the other 
hand, the cobalt e g orbitals interact with a linear combi- 
nation of oxygen orbitals of the e g symmetry, increasing 
the energy of the antibonding a* band. A center of the 
cr*(e s ) band is located by 10D 9 higher than the center of 
the 7r*(t 2 g) band (the same energy splitting as in the CF 
theory) . Note that in the LS state, the a* band is empty. 

The on-site electron interactions as well as distortions 
of the C0O6 octahedron, with the C0O5 pyramid being 
the ultimate case, shift the energy bands in the same 
way as in the case of discrete levels. For IS Co 3+ ions 
in the pyramidal environment, the a* band occupied by 
electrons with the prevailing spin projection is separated 
from the conduction band and shifted down to the top 
of the valence band [as shown in Fig. 30(b)]. On the 
other hand, in the case of the HS state in octahedral 
environment, the occupied a* band is situated higher in 
energy and is twice as large [Fig. 30(c)]. 

Figure 30(b) presents a simplified electronic structure 
originating from a 1:1 mixture of the electronic states of 
C0O6 octahedrons and C0O5 pyramids with the Co 3+ 
ions being in the LS and IS states, respectively. We ex- 
pect this state to be realized in GdBaCo 2 05.5o below the 
metal-insulator transition. The high-temperature elec- 
tronic structure corresponding to HS Co 3+ ions in octa- 
hedral environment and IS ions in pyramidal surround- 
ing is depicted in Fig. 30(c). One can see that the most 
prominent difference between these two pictures is a band 
gap in the first case, while the occupied a* band overlaps 
with the conduction band in the latter one. Whether 
the energy gap emerges or not is determined mainly by 
the position of the a* bands as well as by their widths. 
The electronic structure of GdBaCo 2 05.5o discussed here 
offers a natural, though much simplified, explanation of 
the metal-insulator transition that is observed to coincide 
with the spin-state one at Tmit ~ 360 K. 




FIG. 31: Temperature dependences of the in-plane resistivity 
(a), Hall coefficient (b), effective mobility /z e ff = Rh/p (c), 
and thermoelectric power (d) in GdBaCo2 0s.5o crystals. The 
dashed lines are results of simulation. 



C. Transport properties of GdBaCo2O5.5 

Having obtained a qualitative picture of the electronic 
structure for the parent compound GdBaCo 2 05.so, we 
can now review its transport properties. The tempera- 
ture dependences of the resistivity p{T), Hall coefficient 
-Rh(T), and thermoelectric power S(T) of GdBaCo 2 Os.5o 
are summarized in Fig. 31 together with an "effective" 
mobility of carriers defined as /z e ff = Rn/p- The latter 
would correspond to a true Hall mobility pn, if only one 
type of charge carriers were present in the system. 

Apparently, all the transport properties presented in 
Fig. 31 behave quite coherently At low temperatures, 
both p(T) and Ru(T) grow exponentially upon cooling 
with an activation energy A s=s 70 meV. An activation 
character of the charge transport with the same energy 
gap can be also inferred from the thermoelectric power 
S(T) which is linear in 1/T. On the other hand, when the 
temperature increases above the metal-insulator transi- 
tion at « 360 K, all these physical quantities, p, i?H, and 
S, acquire small and almost temperature-independent 
values: p a \, « p c w 600 /ilJcm, i? H ~ +1.8xl0~ 5 cm 3 /C, 
and S ~ —4 pV/K. This small Hall coefficient would cor- 
respond to an unrealistically high density of carriers of 
w 20 holes per Co ion, if one assumes only one type of car- 
riers to participate in the charge transport; correspond- 
ingly, the apparent mobility /i e ff becomes unrealistically 
small above Tmit [Fig. 31(c)]. One can naturally con- 
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elude, therefore, that the Hall signal is strongly reduced 
at high temperatures because both electrons and holes 
move almost equally and their contributions to i?n, be- 
ing very close in magnitude but opposite in sign, almost 
cancel each other. The metallic charge transport pro- 
vided by both electrons and holes is presumably respon- 
sible for the small value of the Seebeck coefficient as well. 
Note that the sign of S at high temperatures is opposite 
to that of the Hall resistivity, which can easily happen 
when different types of charge carriers are providing del- 
icately balanced contributions: The holes' contribution 
to i?H is a bit larger than the electrons' one, while in S 
the balance appears to be opposite. 

One more interesting feature in Figs. 31(b) and 31(d) 
is that both Ru and S change their sign from nega- 
tive to positive at T w 200 K, long before approaching 
the metal-insulator transition. This sign change indi- 
cates that the electrons and holes in GdBaCo205.so keep 
providing comparable contributions to the conductivity 
down to quite low temperatures, and thus keep compet- 
ing for the dominating role in the Hall response and ther- 
moelectric power. On one hand, this is not surprising, 
given that the neutrality condition in a non-doped semi- 
conductor requires the number of electrons to match the 
number of holes, n = p. On the other hand, however, the 
switch in sign of Rn and S would imply that the relative 
mobilities of electrons and holes should quickly change, 
with their ratio fi n /^p being ;§> 1 at low temperatures 
but Hn/fip <C 1 at high temperatures. 

One possibility is that such asymmetric mobility evo- 
lution actually takes place, though we cannot point to 
any obvious reason for that. Another possibility to be 
considered is that the mobility ratio is essentially tem- 
perature independent, but instead there are more than 
two types of carriers in the system. Namely, there 
may be two distinct types of holes (the total number 
of holes p = pi + P2 = n) with very different mobil- 
ities fj,pi <C /i n <C fJ- P 2 and different activation ener- 
gies. Whenever the temperature is modified, the relative 
shares of the two kinds of holes should change, and their 
average mobility should change accordingly, mimicking a 
gradual mobility evolution. 

In order to understand where the different types of car- 
ries may be coming from, we can refer to the qualitative 
electronic structure of GdBaCo 2 05 50 sketched in Fig. 
30(b). The bands nearest to the Fermi level are the con- 
duction band and a narrow band composed of occupied 
7r* and a* cobalt d-bands, with Ep being located exactly 
in the center of the band gap at T = 0. Apparently, at 
the lowest temperatures, only excitations from the top of 
the d-band to the bottom of the conduction band may 
take place, and thus the carriers in these bands should 
govern the transport properties. It is reasonable to ex- 
pect the mobility of holes in the narrow ci-band composed 
of almost localized states (p p i) to be much lower than 
that of electrons in the conduction band (/i„), though 
the latter is not very high either, since the conduction 
band also includes a* cobalt states (unoccupied). Under 



these conditions, the electrons thermally activated into 
the conduction band should make a predominant contri- 
bution to the low-temperature charge transport, provid- 
ing a simple activation behavior for the resistivity and the 
Hall coefficient, p « enfi n and i?H ~ — 1/en. The ther- 
moelectric power for electrons in the conduction band, 
S n , should follow a conventional law: 
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where the energy Eq corresponds to the conduction-band 
bottom, and r is determined by the scattering mecha- 
nism (through the energy dependence of the relaxation 
time t ~ E r ). Note that a similar expression is valid 
for the thermoelectric power of holes, S p , as well, but 
their low mobility does not allow them to compete with 
electrons at low temperatures, because the total thermo- 
electric power for the two-carrier system is calculated, in 
the simplest model, as 
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where conductivity of electrons, <j n (as well as holes, a p ) 
is directly proportional to their mobility. 

Upon increasing temperature, a wider energy range of 
electronic states becomes available for thermal excita- 
tion, and at some point, a number of electrons excited 
from the valence band to the conduction band should 
become considerable. While the behavior of electrons in 
the conduction band is the same regardless of where they 
are activated from, the holes generated in the valence 
band clearly differ from those in the d band, because the 
valence band is much wider than the d band and the mo- 
bility is expected to be larger in the wider band. As soon 
as the valence-band holes join the transport, the share 
of holes in the conductivity, Hall coefficient, and ther- 
moelectric power quickly increases, and at some temper- 
ature, it exceeds the contribution of electrons owing to 
the high mobility /i P2 > as a result, Rn(T) and S(T) 
change their sign from negative to positive. 

To illustrate the latter picture, in Fig. 31 we show 
results of a simulation (dashed lines) for p(T), Rn(T), 
fj, e ff(T), and S(T) obtained for the three-band model, 
where the mobility ratios were taken as fi Pl : fj, n : fj, P2 = 
0.1 : 1 : 7. One can see that even such very simplified 
model provides a reasonable description for all the mea- 
sured transport properties on the insulator side of the 
metal-insulator transition. We can conclude, therefore, 
that the electronic structure of GdBaCo20s.5o proposed 
in Sec. IIVBI gives a consistent explanation not only for 
the metal-insulator transition, originating from the spin- 
state transition, but also can account for the behavior of 
resistivity, Hall coefficient and thermoelectric power as 
well. 
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D. The origin of giant magnetoresistance 

At this point, it is interesting to consider the mech- 
anism responsible for the giant magnetoresistance ob- 
served in GdBaCo 2 05 +x and other isostructural com- 
pounds. Although the MR behavior in GdBaC^Os+z 
bears some resemblance to that in manganites (in par- 
ticular, the MR in both systems is related to the AF-FM 
transition), there are several dissimilar features pointing 
to different MR mechanisms. Most important is that 
GdBaCo205+a;, in contrast to manganites, remains to 
be an insulator even in the FM state, indicating that 
neither the double-exchange mechanismi 2 ^ nor perco- 
lation through some metallic phase is relevant here^ the 
MR does not seem to originate from the "spin- valve" ef- 
fects, i.e., from a tunneling between two spin-polarized 
metallic regions, either. Our analysis of the resistivity in 
the AF and FM states (Fig. 25) has shown that the ac- 
tivation energy for carriers diminishes considerably upon 
the AF-FM transition, and it is mainly a change in the 
density of carriers, rather than in their mobility, that 
affects the resistivity. 

By examining the crystal and magnetic structure of 
GdBaCo205.5o (Fig. 29), one can easily find out that 
the charge transport, at least along the c axis, should be 
governed by nominally non-magnetic C0O2 layers, which 
go along the ac plane and are composed of cobalt ions 
with octahedral oxygen environment. It appears, there- 
fore, that the charge-carrier doping of these octahedral 
C0O2 planes is to a large extent determined by the rela- 
tive orientation of the magnetic moments in neighboring 
FM ladders. To distinguish this intriguing MR scheme 
from other MR mechanisms, in a recent papenii we have 
coined a term "magnetic field-effect transistor" , based on 
a hypothetic structure where the charge-carrier injection 
into a 2D semiconducting channel is controlled not by 
an electric field, but by a magnetic state of neighboring 
"ligands" . 

One might wonder how the reorientation of very weakly 
coupled FM ladders may affect the charge transport, 
bringing about a giant MR, particularly in the nominally 
non-magnetic C0O2 ac layers. Qualitatively, the MR 
mechanism appears to be quite simple: GdBaCo205.5o is 
a narrow-gap insulator, where the carrier generation goes 
through formation of electron-hole pairs, i.e., through 
formation of inevitably magnetic Co 2+ and Co 4+ states, 
whose energy must depend on the surrounding magnetic 
orderiii In other words, the energy gap that opens at 
the Fermi level due to the crystal-field splitting and on- 
site interactions must be further modified by the spin- 
dependent exchange interactions with neighboring Co 
ions. As we already know, the carriers generated by ther- 
mal excitations or doping always favor the FM coupling 
between ladders (Fig. 28); thus, the FM arrangement of 
ladders should, in turn, also make the carrier generation 
easier. One can conceive this interrelation in a way that 
the energy for generating carriers in the AF state also in- 
cludes a penalty that should be paid upon frustrating the 



relative order of the antiferromagnetically-coupled lad- 
ders. When a magnetic field aligns the FM ladders, this 
penalty is removed, which lowers the barrier for carrier 
generation. One can expect the reduction of the insulat- 
ing gap upon the field-induced AF-FM transition to be 
of the same order as the inter-atomic exchange energy 
J ~ 25 meV, which is indeed observed experimentally 
(Fig. 25). 

Lastly, an interesting point to be emphasized in the 
MR behavior is a strong amplification effect, whereby 
an apparently small energy of the magnetic field 
g/iB^/fcB ~ 1 K becomes capable of changing the carri- 
ers' activation energy by several hundreds K. This gives 
one more evidence for the cooperative nature of the ob- 
served MR. 



E. Evolution of physical properties of 
GdBaCo205+ x with oxygen content 

The main motif governing the behavior of the 
GdBaCo205+a; compound is obviously a competition of 
various kinds of ordering, involving charge, spin, orbital 
and structural degrees of freedom. While for several par- 
ticular compositions (x ks 0,0.5,0.7, etc.), this ordering 
competition is won by certain homogeneous states, the 
majority of the phase diagram appears to be occupied 
by mixtures of various ordered phases (Fig. 28). 

In the end-point x = compound, which corresponds 
to the 50% electron doping, the Co 2+ and Co 3+ ions or- 
der into alternating lines running along the a axis, where 
each cobalt spin is aligned antiferromagnetically to six 
nearest neighbors to form a G-type AF structure!^ This 
means that in GdBaCo20s.oo all the nearest-neighbor 
superexchange interactions, be they between the Co 2+ 
ions or Co 3+ ions, are antiferromagnetic. Under such 
conditions, one would expect the charge doping, which 
converts some Co 2+ ions into Co 3+ ones, to cause no 
considerable effect on the magnetic behavior; this should 
certainly be the case if just the size of AF-ordered spins 
is modified. In reality, however, the oxygen doping im- 
mediately brings in some sort of ferromagnetism. Since 
GdBaCo205+2; retains its insulating behavior all the way 
up to x ~ 0.5, one has to invoke a qualitative change in 
the orbital order to account for the emerging FM be- 
havior. Apparently, the orbital order that underlies the 
AF spin ordering in GdBaCo20s.oo is stable only for the 
particular unidirectional charge order which is realized at 
the 50% electron doping. When the holes are introduced 
upon oxygen doping, they immediately frustrate the pat- 
tern of the charge and orbital order in C0O2 planes. 
In fact, this doping-induced frustration offers a clue to 
the origin of nanoscopic phase separation, whose spec- 
tacular manifestations appear in the magnetic behavior: 
To avoid frustration, the charge-ordered AF matrix ex- 
pel the doped holes, which are then segregated in small 
droplets of a new FM phase, as sketched in Fig. 32; the 
evidence for this phase segregation is found in the mag- 
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FIG. 32: A sketch of the nanoscopic-phase-separated state in 
GdBaCc^Os+z at oxygen concentrations x ~ 0.1 — 0.2. Upon 
oxygen doping, nanoscopic droplets of the orbital-ordered FM 
(OO-FM) phase emerge in the matrix of the charge-ordered 
antiferromagnet (CO-AF). As x increases, these droplets grow 
and should reach percolation at x » 0.25, when a half of the 
AF phase is switched into the FM one. Note that even in 
small droplets, the FM moments should be aligned along one 
of the in-plane tetragonal axis because of the strong spin-orbit 
coupling. 

netization data, as we discussed in Section III D. 

An electronic phase separation, arising from the ten- 
dency of an ordered state to avoid frustration, is quite 
common for strongly correlated systems. Usually, the 
frustration is induced by the charge carrier motion over 
an AF insulating background; the system, therefore, 
tends to segregate the mobile carriers into metallic re- 
gions with a modified magnetic order. as What is pecu- 
liar in GdBaCo205.oo is that the AF and FM phases here 
are both insulating; in this case, they differ only in the 
charge and orbital order. Presumably, these phases cor- 
respond to the charge-ordered insulator (50% electron 
doping) and the undoped band insulator (zero doping). 

An interesting issue is whether the oxygen ions partic- 
ipate in the phase separation or it is mostly an electronic 
phenomenon. Based on the very small size of the FM- 
phase droplets - they are not detected by the X-ray scat- 
tering fSec. HIIA"|) . and exhibit just a superparamagnetic 
behavior at moderate oxygen doping - we can conclude 
that the oxygen ions are not involved in the large-scale re- 
distribution. Though the negatively-charged oxygen ions 
may rearrange locally, being dragged by the Coulomb 
interaction with the doped carriers, this motion is ap- 
parently insufficient to cause the domains to grow. Con- 
sequently, the observed nanoscopic phase separation at 
x < 0.5 is most likely governed by the electronic inter- 
actions in the C0O2 planes and thus may be a generic 
feature of layered cobalt oxides. 

The phase-separated state depicted in Fig. 32 gives 
a consistent explanation for the magnetic and transport 
properties of GdBaCo20s +2: in the electron-doped region 
(0 < x < 0.5). As long as the FM droplets are small and 
located far apart, they exhibit a superparamagnetic be- 
havior with thermomagnetic irreversibility features orig- 
inating from the strong spin-orbit coupling. With in- 
creasing oxygen content, the FM droplets grow until they 
achieve a percolation within the C0O2 planes at x m 0.25; 
from that point, the initial AF phase is kept in inclusions 
embedded in the FM matrix. Eventually, the system be- 
comes homogeneous again as the composition approaches 
that of the undoped parent insulator, x = 0.50. Given 



that both the FM and AF phases are insulators, it is not 
surprising that the conductivity of crystals is virtually 
independent on doping in the entire electron-doped re- 
gion, and that it goes through variable-range hoping of 
localized carriers. 

While in the electron-doped region, GdBaCo205 +:r 
tends to phase separate into two insulating components, 
namely, the undoped and 50%-electron doped phases, the 
hole doping (x > 0.5) also quickly destabilizes the ho- 
mogeneous parent state. However, the phase separation 
in the hole-doped region bears several important differ- 
ences from the behavior of electron-doped crystals. First, 
the scale of resulting domains is much larger, implying 
that the oxygen ions are involved in the macroscopic 
(or mesoscopic) rearrangement process and redistribute 
over the crystal together with the doped carriers; other- 
wise, the long-range Coulomb interactions would prevent 
the formation of such a macroscopically inhomogeneous 
state. Second, the nearest stable phase is located at a 
considerably lower doping level of ~ 18 — 20% of holes 
(x « 0.68 — 0.70). And finally, the most important is that 
the stable hole-doped phase appears to be metallic, in 
contrast to insulating electron-doped compositions. The 
latter explains the remarkable asymmetry in the doping 
dependence of the conductivity in Fig. 7. 

Now, we can only speculate on the nature of the metal- 
lic hole-doped phase, since a detailed information on the 
crystal structure of GdBaCo205+ :c at x « 0.68 — 0.70 is 
still lacking. Apparently, that phase should be well struc- 
turally ordered, since one can hardly imagine so sharp a 
FM-AF transition, as shown in Fig. 18(c), to takes place 
in a system with a significant disorder. Therefore, some 
sort of oxygen ordering is expected to occur in GdO^ 
layers, similar to the formation of perfect oxygen chains 
in the parent compound with x = 0.50; one can imme- 
diately find out that the composition x « 0.68 — 070 
best fits the "Ortho-IE" superstructure where two filled 
oxygen chains in GdOz layers alternate with one empty 
chain. 

If the Ortho-Ill structure is actually realized in 
GdBaCo 2 5+3; with x — 0.7, and the FM-AF transition 
has the same purely magnetic origin as in the parent 
compound with x — 0.50, one inevitably arrives at the 
conclusion that some kind of weakly coupled FM ladders 
or layers are formed in the x ~ 0.7 phase as well. One 
can further speculate that the double octahedral C0O2 ac 
layers are still non-magnetic in this Ortho-Ill phase, and 
the FM ladders are just located farther apart. The main 
difference with the parent x = 0.50 composition is that 
for x ~ 0.7, about 20% of cobalt ions should acquire the 
Co 4+ state. If all the holes are gathered in the FM layers, 
where they can readily move, those layers would possess 
a quarter-filled metallic band. Note that on the macro- 
scopic scale, the conductivity would exhibit a true metal- 
lic behavior only if the crystal has a perfectly ordered 
oxygen subsystem; otherwise, the metallic ladders/layers 
would be always interrupted by non-conducting regions. 
In the latter case, a variable-range hopping between 
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metallic fragments should be expected^ The measured 
resistivity of x = 0.70 crystals indeed exhibits an un- 
usual behavior, tending to saturate at low temperatures 
at finite but surprisingly large values (Fig. 6). 

While the suggested very speculative picture of the 
metallic x ~ 0.7 phase can easily account for the ob- 
served magnetic and transport properties, more work is 
obviously necessary to clarify whether it has anything to 
do with the actual microscopic state. In any case, this 
intriguing phase with self-organized metallic paths bears 
a clear resemblance to the stripe structures in high-T c 
cupratesfii and the phase-separated conducting phases in 
manganites, 2 and thus its study may help in understand- 
ing the unusual transport in other oxides as well. 



V. SUMMARY 

In order to gain insight into the general behavior of 
square-lattice layered cobalt oxides, we have performed a 
systematic study of the transport, magnetic, and thermo- 
electric properties of GdBaCo205 +a; single crystals over 
a wide range of oxygen content, < x < 0.77. The high- 
quality crystals were grown by the floating-zone method 
and their oxygen content was precisely tuned to required 
values by successive annealing treatments so that the en- 
tire phase diagram could be spanned; particularly, in the 
critical regions of the diagram, the composition was mod- 
ified with steps as small as Ax = 0.001 — 0.01. 

The "parent" compound GdBaCo20s.5o - the compo- 
sition with all cobalt ions in the Co 3+ state - is shown to 
be a metal at high temperatures, but it switches into 
a band insulator (a narrow-gap semiconductor) below 
the metal-insulator transition at Tmit ~ 360 K, which 
is manifested in step-like anomalies in all the measured 
properties, p(T), S(T), Ra(T), and M(T). The gap 
opening at the Fermi level is associated with a coopera- 
tive spin-state transition, which seems to involve a half of 
the cobalt ions, while the other half keeps its spin state 
unchanged. Two non-equivalent cobalt positions are cre- 
ated by the oxygen ions which order in GdOo.s layers 
into alternating filled and empty rows running along the 
a axis. Consequently, the C0O2 planes also develop a 
spin-state order consisting of alternating rows of Co 3+ 
ions in the S — 1 and S — states. The overall mag- 
netic structure of GdBaCo20s.5o below Tmit may be con- 
ceived of as a set of magnetic "2-leg ladders" , which are 
composed of S = 1 Co 3+ ions and are separated from 
each other by non-magnetic layers. Upon cooling below 
Tq ~ 300 K, these ladders develop a ferromagnetic order; 
however, whether a macroscopic magnetic moment will 
emerge or not after the FM ladders are formed depends 
on inter-ladder coupling. According to the magnetization 
measurements performed on detwinned single crystals, 
the cobalt spins exhibit a remarkably strong Ising-likc 
anisotropy, being aligned along the oxygen-chain direc- 
tion; thus the FM ladders in GdBaCo205.5o can only 
form a relative ferromagnetic or antiferromagnetic ar- 



rangement. Owing to the very weak strength of this 
coupling between ladder, their relative magnetic order 
can be easily altered by temperature, doping, or mag- 
netic fields, bringing about FM^AF transitions and a 
giant magnetoresistance. 

A study of GdBaCo205+a; crystals with the oxygen 
content deviating from x = 0.50 has revealed that this 
layered compound is a unique filling-control system that 
allows a continuous ambipolar doping of the C0O2 planes, 
that is, the doping level can be continuously driven across 
the parent insulating state. This continuous doping is 
manifested in spectacular singularities of the transport 
properties as the system approaches the undoped state: 
the thermoelectric power, for instance, tends to diverge 
at x — 0.50, and abruptly changes its sign upon crossing 
this peculiar point. 

As soon as the C0O2 planes are doped with more than 
a few percent of charge carriers, be they electrons or 
holes, their homogeneous state is found to become unsta- 
ble, and the system exhibits a very strong tendency to 
separation into several peculiar ordered phases. The re- 
sulting nano- or mesoscopic phase mixture is clearly man- 
ifested in the transport and magnetic properties of the 
crystals. For the entire electron-doped region (x < 0.5), 
GdBaCo205+a; tends to phase separate into two insu- 
lating components, the undoped phase and the charge- 
ordered AF insulator with 50% electron doping. As a re- 
sult, the conduction in GdBaCo205+2; crystals is almost 
doping independent in the entire range < x < 0.45, 
and goes through a variable range hopping of electrons. 
In the hole-doped region, the C0O2 planes also expel the 
doped carriers, protecting the undoped domains, but the 
second phase generated upon the phase separation turns 
out to be metallic, in contrast to the electron-doped crys- 
tals. Consequently, the GdBaCo20s+ x crystals gradually 
evolve towards a metallic state with hole doping, and thus 
the doping dependence of the conductivity exhibits a re- 
markable asymmetry with respect to whether the C0O2 
planes are doped with electrons or holes. 

Although a conclusive picture of the metallic hole- 
doped phase has not been established yet, such metallic 
state self-organized within the electronically inhomoge- 
neous system bears a clear resemblance to the behav- 
ior of CMR manganites and high-T c cuprates, and thus 
its study may offer a clue to understanding the unusual 
charge transport in transition- metal oxides. 

To conclude, the GdBaCo205+ x compound has proved 
to be an interesting system with a very rich phase dia- 
gram originating from the competition of various spin- 
charge-orbital ordered phases. What turns it into a 
model system and makes it particularly interesting for 
a detailed study is the possibility of continuous ambipo- 
lar doping of C0O2 planes and the exceptionally strong 
cobalt-spin anisotropy, which dramatically narrows down 
the range of possible spin arrangements. 
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